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ABSTRACT 

We  present  a  sample  of  386  BL  Lacertae  candidates  identified  from  2860  deg2 8 9  of  the  Sloan  Digital  Sky  Survey 
(SDSS)  spectroscopic  database.  The  candidates  are  primarily  selected  to  have  quasi-featureless  optical  spectra  and 
low  proper  motions  as  measured  from  SDSS  and  USNO-B  positions;  however,  our  ability  to  separate  Galactic  from 
extragalactic  quasi-featureless  objects  (QFOs)  on  the  basis  of  proper  motion  alone  is  limited  by  the  lack  of  reliable 
proper  motion  measurements  for  faint  objects.  Fortunately,  high  proper  motion  QFOs,  mostly  DC  white  dwarfs, 
populate  a  well-defined  region  of  color  space,  approximately  corresponding  to  blackbodies  with  temperatures  in  the 
range  7000-12,000  K.  QFOs  with  measurable  redshifts  or  X-ray  or  radio  counterparts  (i.e.,  evidence  of  an  extra- 
galactic/active  galactic  nucleus  nature)  loosely  follow  a  track  in  color  space  that  corresponds  to  power-law  continua 
plus  host  galaxy  starlight,  with  typical  power-law  slopes  in  the  range  1  <  aopt  <2(f,  <x  v~a).  Based  largely  on  this 
remarkably  clean  color  separation,  we  subdivide  the  sample  into  240  probable  candidates  and  146  additional  less 
probable  (likely  stellar)  candidates.  The  probable  BL  Lac  candidates  have  multiwavelength  properties  consistent 
with  the  range  of  previously  known  BL  Lac  objects,  with  an  apparent  preponderance  of  objects  with  synchrotron 
peaks  at  relatively  high  energies  (HBL/XBL  type).  The  majority  of  the  154  objects  with  measurable  redshifts  have 
z  <  1,  with  a  median  of  0.45;  there  are  also  a  handful  ofhigh-redshift  objects  extending  up  toz  =  5.03.  We  identify  a 
small  number  of  potential  radio-quiet  BL  Lac  candidates,  although  more  sensitive  radio  observations  are  needed  to 
confirm  their  radio-quiet  nature. 
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1.  INTRODUCTION 

Among  the  rarest  and  most  extreme  of  the  various  observa¬ 
tional  classes  of  active  galactic  nuclei  (AGNs)  are  the  BL  Lacertae 
objects.  The  defining  traits  of  BL  Lac  objects  generally  include 
a  broad,  rapidly  variable,  nonthermal  spectral  energy  distribution 
(SED)  with  significant  power  from  radio  wavelengths  to  X-rays 
and  y-rays;  the  lack  of  the  strong  optical  emission  lines  that  typify 
other  AGNs;  high  and  variable  polarization  at  optical  and  radio 
wavelengths;  and  a  compact,  core-dominated  radio  morphology 
(Kollgaard  1994).  All  of  the  above  properties  can  be  reproduced 
naturally  by  a  model  in  which  the  dominant  source  of  the  ob¬ 
served  emission  is  a  relativistic  jet  of  plasma  directed  within  a 
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small  angle  to  the  line  of  sight  (e.g.,  Blandford  &  Rees  1978). 
This  picture  is  particularly  appealing  because  there  is  a  known 
class  of  object  that  plausibly  consists  of  physically  identical  sys¬ 
tems  in  which  the  jet  is  directed  at  a  more  substantial  angle  to  the 
line  of  sight,  namely,  the  Fanaroff-Riley  class  I  (FR I;  low  power; 
Fanaroff  &  Riley  1974)  radio  galaxies.  The  beamed  counterparts 
of  more  powerful  FR  II  radio  galaxies  are  generally  believed  to 
be  flat-spectrum  radio  quasars  (FSRQs),  although  FSRQs  have 
emission-line  properties  more  like  those  of  typical  broad-line 
AGNs.  BL  Lac  objects  (at  lower  typical  luminosity)  and  FSRQs 
(at  higher  luminosity)  are  often  grouped  together  under  the  desig¬ 
nation  of  blazars.  This  scenario  is  favored  in  unification  schemes 
for  radio-loud  AGNs  (Urry  &  Padovani  1995). 

Traditionally,  BL  Lac  objects  have  been  discovered  from  ob¬ 
servations  at  either  radio  (e.g.,  Shekel  et  al.  1991)  or  X-ray  (e.g., 
Stocke  et  al.  1991;  Schachter  et  al.  1993;  Perlman  et  al.  1996) 
frequencies.  Such  radio-  and  X-ray-selected  samples  of  BL  Lac 
objects  are  known  to  be  composed  of  objects  with  systematically 
different  properties,  such  as  the  peak  frequencies  of  the  SEDs,  as 
discussed  by  Padovani  &  Giommi  (1995).  Those  authors  coined 
the  terms  “low-energy  cutoff  BL  Lac  objects”  (LBLs11)  and 
“high-energy  cutoff  BL  Lac  objects”  (HBLs)  to  distinguish  the 
types  of  objects  predominantly  found  in  radio  and  X-ray  se¬ 
lected  samples,  respectively;  the  obvious  acronyms  RBL  and 
XBL  are  often  used  to  make  the  same  distinction.  It  has  been 
proposed  that  blazar  properties  follow  a  continuous  sequence 
from  BL  Lac  objects  into  the  regime  of  FSRQs;  in  this  scenario 


11  In  more  recent  context  the  word  “cutoff,”  which  refers  to  a  break  in  the 
energy  distribution  of  the  relativistic  electrons  responsible  for  the  synchrotron 
emission,  is  often  replaced  with  “peaked,”  taken  to  refer  to  the  energy  corre¬ 
sponding  to  the  synchrotron  peak. 
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an  object’s  observed  characteristics  are  determined  by  the  angle 
between  the  jet  axis  and  the  line  of  sight,  the  overall  jet  power, 
and  possibly  another  parameter  (Fossati  et  al.  1998;  Ghisellini 
et  al.  1998). 

Large  optical  surveys,  such  as  the  Sloan  Digital  Sky  Survey 
(SDSS;  York  et  al.  2000)  and  the  Two-Degree  Field  QSO  Red- 
shift  Survey  (2QZ;  Boyle  et  al.  2000),  have  the  potential  to  reveal 
new  populations  ofBL  Lac  objects  with  different  properties  from 
the  well-known  radio-  or  X-ray-selected  objects.  Londish  et  al. 
(2002)  presented  a  sample  of  several  dozen  objects  from  the  2QZ 
selected  on  the  basis  of  having  quasi-featureless  optical  spectra 
resembling  those  of  BL  Lac  objects;  however,  some  of  these 
objects  were  found  to  be  significantly  fainter  in  the  radio,  relative 
to  the  optical,  than  classical  BL  Lac  objects.  At  least  one  of  these 
latter  objects  has  been  confirmed  to  be  extragalactic  via  a  redshift 
measurement  obtained  in  a  follow-up  observation  (Londish  et  al. 
2004);  along  with  a  small  number  from  other  searches  (Fan  et  al. 
1999;  Anderson  et  al.  2001),  these  objects  are  potential  radio¬ 
quiet  BL  Lac  objects  or  weak-lined  radio-quiet  quasars.  Whether 
or  how  such  objects  fit  into  the  blazar  sequence  or  more  generally 
into  the  AGN  unification  picture  are  open  questions. 

In  the  present  work,  we  describe  a  sample  of  386  optically 
identified  BL  Lac  candidates  compiled  from  2860  deg2  of  the 
SDSS  spectroscopic  survey.  As  in  Londish  et  al.  (2002),  these 
objects  are  selected  to  have  quasi-featureless  optical  spectra  and 
are  identified  solely  on  the  basis  of  their  optical  properties  and 
lack  of  significant  proper  motions.  We  discuss  the  degree  to  which 
the  sample  remains  contaminated  by  featureless  stars,  those  ei¬ 
ther  lacking  reliable  proper  motion  measurements  or  simply  hav¬ 
ing  proper  motions  below  the  detection  threshold.  We  subdivide 
the  sample  into  240  “probable”  and  146  “possible”  candidates 
based  on  all  available  information.  We  search  the  sample  for  ex¬ 
ceptional  objects,  such  as  radio-quiet  BL  Lac  objects  (or  lineless 
radio-quiet  quasars),  and  explore  how  our  sample  fits  into  the 
overall  understanding  of  the  blazar  class.  Our  probable  BL  Lac 
candidates  are  not  a  complete  sample  but  are  more  numerous 
than  any  previous  sample  of  BL  Lac  objects  of  which  we  are 
aware  and  represent  a  significant  addition  to  the  ^400  previ¬ 
ously  known  BL  Lac  objects. 

In  §  2  we  describe  the  SDSS  observations  and  data  processing. 
In  §  3  we  discuss  the  selection  of  quasi-featureless  spectra  and  the 
use  of  proper  motion  to  reject  weak- featured  stellar  objects.  In  §  4 
we  cross-correlate  our  sample  with  the  ROSAT  All-Sky  Survey 
(RASS)  X-ray  source  catalogs  and  with  the  Faint  Images  of  the 
Radio  Sky  at  Twenty  cm  (FIRST)  and  NRAO  VLA  Sky  Survey 
(NVSS)  radio  catalogs.  In  §  5  we  analyze  the  redshift  distribu¬ 
tion,  optical  colors,  and  multiwavelength  properties  of  our  BL 
Lac  candidates.  In  §  6  we  lay  out  some  prospects  for  follow-up 
studies,  and  in  §  7  we  summarize  our  primary  results. 

2.  SDSS  OBSERVATIONS  AND  DATA  PIPELINE 

The  SDSS  is  a  multi-institutional  effort  to  obtain  high- 
quality  five-band  (SDSS  ugriz;  Fukugita  et  al.  1996;  Hogg  et  al. 
2001;  Smith  et  al.  2002)  photometric  imaging  of  a  large  fraction 
(« 10,000  deg2)  of  the  sky,  principally  the  northern  Galactic  cap 
and  a  stripe  along  the  celestial  equator,  with  spectroscopic  follow¬ 
up  of  interesting  samples  of  objects,  primarily  galaxies  (~106) 
and  quasars  (~105).  The  SDSS  uses  a  dedicated  2.5  m  telescope 
at  Apache  Point  Observatory  with  a  2?5  30  CCD  mosaic  camera 
(Gunn  et  al.  1998)  and  two  multiobject  spectrographs  with  320  3" 
fibers  each  (e.g.,  Stoughton  et  al.  2002).  Imaging  data  are 
processed  through  automated  software  pipelines  ( Lupton  et  al. 
200 1 ;  Pier  et  al.  2003)  in  which  the  data  are  photometrically  and 
astrometrically  calibrated  and  objects  are  detected,  measured. 


and  classified.  The  photometric  calibration  is  typically  better 
than  2%  (Ivezic  et  al.  2004),  and  the  astrometric  calibration  is 
better  than  0"1  rms  per  coordinate. 

2. 1 .  Targeting  for  Spectroscopy 

Once  a  region  of  sky  has  been  imaged,  subsets  of  the  objects 
detected  are  then  selected  for  spectroscopic  observations  on  the 
basis  of  their  photometric  and  other  properties.  Resolved  objects 
are  targeted  to  r  =  17.77  for  the  main  galaxy  sample  (Strauss 
et  al.  2002),  and  luminous  red  galaxies  are  targeted  based  on  their 
colors  to  r  =  19.5  (Eisenstein  et  al.  2001).  Quasar  candidates  are 
targeted  as  outliers  from  the  stellar  locus  or  as  pointlike  FIRST 
sources  to  i  =  19.1  for  the  low-redshift  sample;  objects  very  red 
in  u  —  g,  g  —  r,  or  r  —  i  (and  having  colors  consistent  with  Lya 
absorption  cutting  into  the  observed  frame  optical  flux)  are  tar¬ 
geted  to  i  =  20.2  for  the  high-redshift  quasar  sample  (Richards 
et  al.  2002).  Ultraviolet  excess  (UVX;  dereddened  u  —  g  <  0.6) 
objects  with  nonstellar  colors  are  explicitly  included  in  the  low- 
redshift  sample.  Objects  can  also  be  targeted  by  a  variety  of 
serendipity  algorithms  to  g  =  20.5  or  i  =  20.5.  These  selection 
routines  cover  regions  of  color  space  outside  the  main  stellar 
locus  but  not  included  by  quasar  target  selection;  they  also 
include  objects  with  radio  counterparts  detected  by  the  FIRST 
survey  (Becker  et  al.  1995;  Ivezic  et  al.  2002)  or  X-ray  counter¬ 
parts  from  the  RASS  ( Voges  et  al.  1999;  Anderson  et  al.  2003). 
Spectroscopic  fibers  are  allocated  by  a  tiling  algorithm  ( Blanton 
et  al.  2003),  which  gives  top  priority  to  the  main  quasar  and  galaxy 
samples  (along  with  fibers  needed  for  calibration  purposes). 

2.2.  Processing  of  Spectra 

Raw  spectral  data  are  reduced  and  calibrated  (e.g.,  Stoughton 
et  al.  2002),  and  individual  spectra,  which  cover  the  wavelength 
range  3800-9200  A  with  an  approximate  resolution  of  2/AA  ss 
1800,  are  then  automatically  processed  to  classify  them  (star, 
galaxy,  or  quasar)  and  determine  redshifts. 

In  our  analysis,  we  use  the  output  of  the  classification/redshift- 
finding  routine  “SpecBS”  (written  by  D.  Schlegel).  This  algo¬ 
rithm  classifies  spectra  by  performing  separate  x2-fits  to  the 
spectra  using  linear  combinations  of  different  models — sets  of 
galaxy  or  quasar  eigenspectra  or  stellar  templates — plus  second- 
or  third-order  polynomials,  with  redshift/radial  velocity  as  a  free 
parameter.  An  object  is  classified  as  a  star,  galaxy,  or  quasar 
according  to  which  type  of  model  produces  the  best  ff  value. 
Redshift  warning  flags  are  set  if  there  are  multiple  models  that 
provide  nearly  equivalent  \2  values  or  if  other  problems  are  pres¬ 
ent,  such  as  limited  wavelength  coverage  due  to  lost  data,  un¬ 
reasonable  behavior  like  negative  templates,  or  large  numbers 
of  statistically  anomalous  data  points. 

The  polynomial  component,  which  may  be  positive  or  neg¬ 
ative,  is  included  to  account  for  (generally  low  level)  variations 
of  individual  objects  from  the  templates  and  to  compensate  for 
contamination  (e.g.,  star-galaxy  blends)  or  flux  calibration  prob¬ 
lems  (cf.  discussion  in  Abazajian  et  al.  2004),  as  well  as  to  assist 
in  fitting  the  spectra  of  unusual  (e.g.,  heavily  reddened)  objects. 
Objects  with  the  quasi-featureless  spectra  we  seek  generally  re¬ 
quire  strong  polynomial  components  in  order  to  achieve  ade¬ 
quate  fits,  as  discussed  in  §  3.  In  the  most  extreme  cases,  the 
polynomial  components  can  completely  dominate  the  fits;  for 
such  objects  the  spectral  classifications  and  redshifts  are  clearly 
meaningless. 

3.  SAMPLE  SELECTION 

In  order  to  avoid  the  selection  biases  that  enter  into  the  con¬ 
struction  of  radio-  and  X-ray-selected  samples  of  BL  Lac  objects 
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to  the  greatest  degree  possible,  we  confine  our  selection  criteria 
strictly  to  the  optical  properties  of  objects  in  the  SDSS  spec¬ 
troscopic  database.  Of  course,  objects  enter  into  that  data  set  via 
a  variety  of  target  selection  algorithms  as  described  in  §  2,  some 
of  which  are  based  on  X-ray  or  radio  emission.  However,  by 
keeping  track  of  the  reasons  for  which  each  object  was  targeted 
for  spectroscopy,  we  can  account  after  the  fact  for  selection 
biases  introduced  by  these  factors  (see  discussion  in  §  5.5). 

Most  of  the  sky  coverage  of  SDSS  is  imaged  in  one  epoch  only. 
Spectroscopic  observations  of  a  given  field  typically  lag  behind 
imaging  by  anywhere  from  approximately  1  month  to  1  yr.  A 
detailed  study  of  the  variability  of  spectroscopically  confirmed 
quasars  using  the  photometric  +  spectroscopic  observations  is 
presented  by  Vanden  Berk  et  al.  (2004).  For  our  purposes,  how¬ 
ever,  (usually)  having  only  two  epochs  of  observation  in  total 
for  a  given  spectroscopic  object  precludes  a  complete  selection 
on  the  basis  of  variability.  Another  optical  property  that  could 
be  used  to  identify  BL  Lac  objects  would  be  polarization  (e.g., 
Jannuzi  et  al.  1993),  but  this  is  not  among  the  capabilities  of  the 
SDSS. 

Therefore,  the  spectroscopic  phase  of  our  BL  Lac  candidate 
selection  is  primarily  based  on  a  single  property:  lack  of  strong 
features  (i.e.,  emission  or  absorption  lines)  in  the  optical  spectra. 
While  this  approach  is  clearly  biased  toward  objects  that  are 
dominated  by  the  continuous  emission  from  the  jet  at  optical 
wavelengths  (and  subject  to  the  effects  of  SDSS  spectroscopic 
target  selection),  this  bias  should  differ  from  those  that  affect 
classical  radio  or  X-ray  selected  samples  and  therefore  may  be 
sensitive  to  new  populations  of  BL  Lac  objects  missed  by  pre¬ 
vious  surveys.  This  is  the  same  approach  used  by  Londish  et  al. 
(2002). 

The  parent  sample  for  our  search  is  the  SDSS  spectroscopic 
database  as  of  2002  mid-July,  comprising  over  345,000  indi¬ 
vidual  spectra  and  covering  about  2860  deg2  of  sky  (roughly  the 
area  covered  by  Data  Release  Two  [DR2];  Abazajian  et  al.  2004). 
The  initial  selection  was  performed  using  a  somewhat  older  ver¬ 
sion  of  the  spectroscopic  data  reduction  pipeline,  approximately 
corresponding  to  Data  Release  One  ( DR1 ;  Abazajian  et  al.  2003). 
In  that  version  there  were  (known)  spectrophotometric  problems, 
especially  in  the  3800-4400  A  range,  such  as  artificial  structure 
resembling  emission  or  absorption  features;  we  employed  cau¬ 
tion  in  our  search  so  as  not  to  mistake  these  and  other  data  re¬ 
duction  artifacts  for  real  spectral  structure.  These  problems  have 
been  corrected  to  a  large  degree  in  DR2  and  later  reductions  (cf. 
Abazajian  et  al.  2004).  As  previously  stated,  we  use  the  SpecBS 
spectroscopic  classifications  and  spectral  fits  (contemporaneous 
with  DR2);  all  other  data  presented  in  this  paper  (e.g.,  photom¬ 
etry,  astrometry,  etc.)  are  from  the  official  DR2  reductions,  unless 
stated  to  the  contrary. 

3.1.  Selection  of  Quasi-F ecitureless  Spectra 

The  goal  of  our  search  is  to  extract  all  high-quality  SDSS  spec¬ 
tra  that  are  consistent  with  the  featureless  or  quasi-featureless 
spectrum  of  a  BL  Lac  object,  without  regard  to  the  continuum 
slope.  Depending  on  intrinsic  luminosity,  redshift,  and  the  details 
of  AGN  physics,  these  objects  may  show  (1)  host  galaxy  sig¬ 
natures  such  as  stellar  absorption  features  (e.g.,  Ca  n  H  and  K)  or 
emission  lines  from  associated  H  n  regions  (in  which  case  there 
is  a  good  chance  that  the  image  of  the  object  will  also  show 
spatially  resolved  structure),  (2)  weak  broad  or  narrow  emission 
lines  associated  with  AGN  activity,  (3)  metal  absorption  lines  or 
Lya  forest  absorption  from  intervening  material  along  the  line 
of  sight,  or  (4)  no  features  whatsoever,  to  the  limit  of  the  SDSS 
spectral  sensitivity  (typically  better  than  a  few  angstroms  equiv¬ 


alent  width  for  narrow  features,  for  the  objects  that  satisfy  our 
signal-to-noise  ratio  [S/N]  criteria;  see  §  3.1.1).  Especially  for 
the  truly  featureless  objects,  we  reiterate  that  the  spectral  clas¬ 
sifications  and  redshifts  obtained  from  the  spectroscopic  pipe¬ 
line  will  be  meaningless. 

A  complication  inherent  in  this  approach  is  that  BL  Lac  ob¬ 
jects  are  not  the  only  type  of  objects  known  to  have  featureless 
optical  spectra.  In  particular,  white  dwarfs  ( WDs)  of  class  DC 
share  this  property  (e.g.,  Wesemael  et  al.  1993;  Londish  et  al. 
2002).  This  is  not  an  issue  for  objects  that  fall  into  classes  ( 1 )— (3) 
above.  To  distinguish  true  BL  Lac  objects  of  class  (4)  from  DC 
WDs  without  resorting  to  multiwavelength  properties,  the  best 
available  test  is  to  use  proper  motion  (see  §  3.2  below).  In  what 
follows,  we  refer  to  objects  that  satisfy  our  spectral  selection 
criteria  as  quasi-featureless  objects  (QFOs);  we  shall  reserve  the 
term  “BL  Lac  candidate”  for  objects  that  satisfy  the  additional 
proper-motion  requirement  described  in  §  3.2. 

3.1.1.  Signal-to-Noise  Ratio  Requirement 

In  order  to  restrict  our  analysis  to  objects  for  which  reliable  spec¬ 
tral  classifications  are  possible,  we  first  impose  an  S/N  require¬ 
ment.  We  analyze  the  spectral  flux  density  (ergs  cm~2  s_1  A-1)  in 
three  500  A  wide  regions  centered  at  4750,  6250,  and  7750  A 
in  the  observed  frame;  these  are  chosen  to  be  regions  in  which 
we  find  SDSS  spectra  to  be  of  generally  high  quality,  with  min¬ 
imal  sky  contamination  and  calibration  problems,  and  that  fall 
within  the  g,  r,  and  i  filters,  respectively.  We  define  the  signal  S 
in  each  of  these  regions  to  be  5=  J2i(fxj/a^)/Yli(l  /a'f),  where 
fxj  and  a i  are  the  flux  density  and  estimated  uncertainty  in 
spectral  element  i.  The  error  N  on  this  quantity  is  simply  N  = 
(£;l/ofr1/2;  therefore,  the  S/N  for  each  region  is  S/N  = 
(£!yi,!7o'2)/(£/T/cr2)1/2.  We  keep  only  objects  for  which  S/N  > 
1 00  in  one  or  more  of  these  spectral  regions.  This  corresponds  to 
fiber  magnitudes  of  roughly  g  <  20.5,  r  <  20.3,  or  i  <  19.6  for 
a  typical  spectroscopic  plate. 

3.1.2.  Objects  Classified  as  Stars  or  Galaxies 

Stars  and  galaxies  together  account  for  the  majority  (89%) 
of  all  SDSS  spectroscopic  objects.  Most  of  these  are  well  fitted 
by  stellar  templates  or  combinations  of  galaxy  eigenspectra, 
respectively,  without  the  need  for  significant  polynomial  com¬ 
ponents  in  the  fits  (see  §  2.2).  This  is  true  only  for  objects  that 
have  characteristic  stellar  or  galactic  SEDs  and  absorption  fea¬ 
tures  of  typical  strength;  such  objects  clearly  do  not  qualify  as 
quasi-featureless.  In  order  to  separate  such  objects  from  po¬ 
tential  QFOs,  we  examine  the  relative  importance  of  the  poly¬ 
nomial  component  in  the  fit  in  two  spectral  regions:  3800-4200 
and  8800-9200  A  (in  the  observed  frame;  near  the  blue  and  red 
ends,  respectively,  of  the  SDSS  spectral  coverage).  We  define 
the  ratio  R  =  P/S,  where  P  is  the  average  value  of  the  poly¬ 
nomial  in  a  region  and  S  is  the  spectral  signal  as  defined  in 
§  3.1.1.  IfS/N  <  25  in  both  regions,  the  values  of  Rare  not  very 
meaningful  and  the  test  cannot  be  applied;  such  objects  are 
rejected  from  further  consideration.  (This  is  similar  to  the  S/N 
requirement  described  in  §  3.1.1,  but  much  less  stringent.)  We 
reject  objects  classified  by  SpecBS  as  stars  or  galaxies  for  which 
Rbiue  <  0-6  and  Rred  <  0.5;  the  blue  cutoff  is  more  stringent 
because  of  the  large  number  of  absorption  features  in  this  wave¬ 
length  region  for  many  stellar  types  and  the  presence  of  the 
Ca  ii  H  and  K  break  in  galaxy  spectra.  This  step  eliminates  ap¬ 
proximately  80%  of  stars  and  96%  of  galaxies. 

We  do  not  reject  objects  classified  as  quasars  based  on  the 
values  ofR.  Because  of  the  range  of  quasar  emission  line  strengths 
and  the  power-law  nature  of  quasar  continuum  emission,  some 
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Fig.  1. — Example  weak-featured  SDSS  spectra.  The  smooth  solid  curves 
show  the  polynomial  components  from  the  pipeline  spectral  fits.  The  feature  at 
5577  A,  most  prominent  in  (b)  and  (c),  is  a  residual  from  sky  subtraction.  Two  of 
the  three  objects  were  misclassified  by  the  spectral  pipeline  and  assigned  er¬ 
roneous  redshifts;  in  all  three  cases  redshift  warning  flags  were  set.  The  dotted 
vertical  lines  mark  the  positions  of  apparent  absorption  features,  (a)  No  features 
are  detected;  this  object  survives  the  stellar  absorption  line  cut  described  in 
§  3.1.2  and  enters  the  sample  of  QFOs.  ( b )  He  i  XX  4026,  4471,  and  5876  are 
detected  in  absorption;  thus,  this  object  is  probably  a  DB  WD.  (c)  Shallow  Swan 
bands  of  the  C2  molecule  are  detected  (dotted  lines  mark  their  approximate 
edges);  this  object  is  a  DQ  WD.  Objects  ( b )  and  (c)  are  rejected  as  QFOs. 

quasi-featureless  objects  may  be  well  fitted  by  combinations  of 
quasar  eigenspectra  (through  the  observed  wavelength  range) 
without  the  need  for  a  strong  polynomial  component. 

Objects  classified  as  stars  that  survive  the  previous  cut  are 
then  subjected  to  a  test  for  absorption  features  appropriate  to  the 
SpecBS-fit  stellar  type.  We  test  stars  classified  as  O,  B,  or  A  for 
absorption  lines  from  He  n,  He  i,  and  the  H  Balmer  series.  We 
test  F  and  G  stars  for  the  H  Balmer  series  and  lines  from  Ca  n, 
Mg  1,  and  Nai.  We  test  K  stars  forCai,  Can,  Mg  i,  and  Nai  lines. 
We  test  M  stars  for  Ca  n  and  Na  i  lines  and  TiO  molecular 
absorption  bands.  We  measure  the  significance  of  all  absorption 
features  by  comparing  the  spectral  flux  in  predefined  continuum 
regions  on  either  side  of  the  absorption  to  the  spectral  flux 
within  the  absorption  feature.  Objects  are  rejected  if  two  or 
more  absorption  features  are  detected  with  >3  a  significance; 
we  require  the  detection  of  at  least  two  features  to  decrease  the 
chance  of  inadvertently  rejecting  extragalactic  objects  showing 
absorption  lines  from  intervening  material.  Examples  of  objects 
with  weak  stellar  features  or  featureless  spectra  are  displayed 
in  Figure  1. 

We  test  objects  classified  as  galaxies  for  the  strength  of  the 
Ca  ii  H  and  K  break  (at  the  pipeline  redshift)  if  it  is  present  in 
the  spectra  and  if  no  redshift  warning  flags  have  been  set.  The 
Ca  ii  H  and  K  break  strength  is  defined  as 

C=  0.14  +  0.86(4+  -fx,-)/A+,  (1) 

where/)  _  and //+  are  the  average  flux  densities  in  the  rest- frame 
wavelength  regions  3750-3950  and  4050-4250  A,  respectively 
(e.g.,  Landt  et  al.  2002).  (The  form  of  this  equation  derives  from 
the  fact  that  it  is  historically  measured  in fv  as  in  Bruzual  A.  1 983, 
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Fig.  2. — Example  SDSS  spectra  of  QFOs  with  noticeable  host  galaxy  fea¬ 
tures  of  varying  strengths  and  secure  redshift  measurements.  The  format  is  the 
same  as  in  Fig.  1.  The  dotted  lines  mark  the  positions  of  stellar  absorption  lines 
at  the  appropriate  redshifts.  Objects  (a)  and  ( b )  are  unambiguous  BL  Lac  can¬ 
didates.  Object  ( c )  represents  a  marginal  case  for  a  BL  Lac  object:  it  has  a  Ca  n  H 
and  K  break  strength  of  0.36,  larger  than  the  classical  cutoff  of  0.25  but  within 
our  acceptable  range  of  <0.4.  All  three  objects  enter  the  sample  of  QFOs. 

rather  than  /).)  We  reject  objects  for  which  C  >  0.4  (as  in  Marcha 
et  al.  1996)  and  \C/bC\  >  3,  where  6C  is  the  uncertainty  on  C. 
Stocke  et  al.  (1991)  adopted  a  more  conservative  cut  of  0.25, 
but  we  err  on  the  liberal  side  in  order  not  to  reject  any  intrinsically 
red  BL  Lac  objects  with  significant  host  galaxy  contamination. 
Several  examples  of  objects  with  a  range  of  host  galaxy  con¬ 
tributions  that  survive  the  polynomial  and  Ca  n  H  and  K  cuts  are 
given  in  Figure  2. 

3.1.3.  Objects  Classified  as  Quasars  and  Other  Emission-Line  Objects 

For  all  objects  classified  as  quasars  and  for  objects  classified  as 
stars  or  galaxies  that  survive  the  cuts  described  in  §  3. 1 .2,  we  test 
for  the  presence  of  emission  lines  at  the  pipeline  redshift.  We 
measure  the  equivalent  width  (EW)  of  an  emission  line  by 
comparing  the  integrated  flux  in  a  predefined  wavelength  region 
around  the  center  of  the  line  with  a  linear  approximation  to  the 
continuum  in  that  region  of  the  spectrum.  The  continuum  level  is 
determined  by  examining  the  flux  in  predefined  regions  on  either 
side  of  the  line.  Special  care  is  taken  to  examine  potential  con¬ 
tinuum  regions  for  the  presence  of  absorption  lines  (e.g.,  metal 
lines  from  intervening  material),  which  could  artificially  lower 
the  continuum  level  and  thereby  exaggerate  EW  measurements. 
We  do  not  attempt  to  automatically  measure  the  EW  of  the  Lya 
emission  line  because  of  the  frequently  complex  absorption  to¬ 
ward  shorter  wavelengths.  This  is  to  protect  against  inadver¬ 
tently  rejecting  any  BL  Lac  objects  at  high  enough  redshift  to 
show  absorption  from  the  Lya  forest  (e.g.,  Fan  et  al.  1999). 

We  reject  all  objects  in  which  any  identified  emission  line  has 
a  rest-frame  equivalent  width  (RFEW)  greater  than  5  A  at  1  <j 
(i.e.,  the  1  er  error  on  RFEW  does  not  overlap  5  A).  This  criterion 
is  in  accordance  with  the  classical  definition  of  a  BL  Lac  object 
but  is  rather  strict  by  modem  standards  (e.g.,  Marcha  et  al.  1996) 
in  that  it  may  cause  us  to  reject  objects  with  featureless  AGN 
continua  plus  emission  lines  from  host  galaxy  star  formation. 
However,  we  employ  it  in  order  not  to  be  overwhelmed  by 
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Fig.  3. — Example  SDSS  marginal  emission-line  objects,  similar  to  Figs.  1 
and  2.  Dotted  lines  again  mark  absorption  features,  and  dashed  lines  mark 
emission  features.  All  three  spectra  show  significant  sky  residuals  at  the  red 
ends;  many  of  the  apparent  emission  features  longward  of  about  7500  A  cannot 
be  trusted,  (a)  The  broad  emission-line-like  feature  observed  near  5630  A  may 
be  Mg  ii  22800;  however,  the  identification  is  not  secure,  primarily  because  of 
the  lack  of  other  features,  so  this  object  survives  the  emission-line  cut  described 
in  §  3.1.3  and  enters  the  sample  of  QFOs.  ( b )  The  broad  lines  appear  weak  and/or 
absorbed  in  this  object;  the  C  iv  21549  line  (observed  near  6800  A  and  possibly 
blueshifted  with  respect  to  Lya )  has  RFEW  close  to  5  A  but  with  substantial  uncer¬ 
tainty.  Therefore,  object  ( b )  survives  as  a  QFO.  (c)  O  n  23727  has  RFEW  >  5  A. 
Object  (c)  is  rejected  as  a  QFO. 

galaxies  with  moderate  star  formation  rates,  which  exist  in  the 
SDSS  spectroscopic  data  in  large  numbers. 

A  line  identification  is  considered  secure  only  if  other  spec¬ 
tral  features  are  present  at  the  same  redshift;  thus,  a  line  iden¬ 
tification  implies  that  a  redshift  can  be  determined.  If  a  single 
emission  feature  appears  to  be  present  but  a  line  identification 
cannot  be  unambiguously  determined,  we  do  not  reject  the  object. 
These  criteria  are  designed  so  that  QFOs  (possibly  having  bumpy 
spectra,  e.g.,  due  to  intervening  absorption)  will  not  be  rejected 
no  matter  what  redshift  the  pipeline  assigns  them.  Example  spec¬ 
tra  of  borderline  emission-line  objects  are  shown  in  Figure  3. 

3.1.4.  Manual  Examination  and  Frequent  Contaminants 

The  selection  procedure  described  above  is  performed  by  an 
automated  algorithm;  however,  the  algorithm  is  designed  to  err 
on  the  side  of  caution  so  as  not  to  inadvertently  reject  unusual 
and  interesting  objects.  The  result  is  that  3527  spectra  (1.02%  of 
the  parent  sample)  survive  through  the  automated  stage  and  must 
be  manually  examined.  At  this  stage  we  stringently  reject  re¬ 
maining  emission-line  objects  for  which  unambiguous  line  iden¬ 
tifications  can  be  made  (with  RFEW  >  5A),  objects  with  visibly 
identifiable  stellar  features  near  zero  redshift,  and  objects  with 
data  reduction  problems  or  various  forms  of  contamination 
(e.g.,  star-galaxy  blends).  All  objects  in  Figures  1-3  survived 
until  at  least  the  manual  selection  stage  because  of  the  conser¬ 
vative  design  of  the  automated  algorithm. 

In  addition  to  the  quasi-featureless  objects  we  seek,  several 
distinct  types  of  objects  survive  in  significant  quantities  until  the 
manual  examination  stage.  Among  these  are  quasars  with  rela¬ 
tively  weak  emission  lines  such  as  the  object  shown  in  Figure  3b. 
Broad  absorption  line  (BAL)  quasars  (e.g.,  Flail  et  al.  2002; 


Reichard  et  al.  2003)  also  survive  in  large  numbers.  Their  fre¬ 
quently  complicated  spectra  render  it  difficult  to  define  continuum 
regions  for  measuring  emission-line  EWs;  we  reject  objects  in 
which  BAL  troughs  prohibit  emission-line  measurements.  An¬ 
other  class  of  frequent  contaminants  are  blue,  low  surface  bright¬ 
ness  galaxies.  The  spectra  of  these  objects  often  show  weak  or 
washed-out  absorption  features;  however,  they  are  easily  identi¬ 
fied  by  their  image  properties.  A  final  type  offrequent  contaminant 
are  white  dwarfs  with  weak-to-undetectable  absorption  features; 
many  of  these  objects  survive  through  the  manual  examination 
step  and  are  discussed  in  the  following  subsection.  Among  weak- 
featured  WDs  several  types  stand  out,  including  weak  DQ  (car¬ 
bon)  WDs  (e.g.,  Liebert  et  al.  2003)  and  weak  DB  (helium)  WDs. 
Kleinman  et  al.  (2004)  present  a  catalog  of  WDs  from  SDSS, 
including  many  with  very  weak  spectral  features  or  apparently 
featureless  spectra. 

3.2.  Removing  Contamination  by  Weak-featured  White  Dwarfs 

After  the  final  examination  of  all  remaining  spectra,  we  are  left 
with  a  sample  of  580  QFOs  out  of  the  original  ~345,000.  Of 
these,  many  are  expected  to  be  DC  WDs,  which  also  fulfill  the 
requirement  of  having  featureless  optical  spectra.  (A  total  of  124 
of  the  580  QFOs  appear  in  the  Kleinman  et  al.  2004  catalog;  cf. 
§§  3.2.3  and  5.3.)  In  addition,  there  may  be  other  weak- lined  stel¬ 
lar  objects  and  possibly  extragalactic  non-BL  Lac  QFOs  (e.g.,  Fan 
et  al.  1999;  Londish  et  al.  2004),  which  would  be  of  great  interest. 

In  order  to  separate  Galactic  from  extragalactic  objects,  several 
tests  can  be  performed.  One  is  the  determination  of  a  redshift/ 
radial  velocity  based  on  spectral  features.  A  second  property 
that  can  indicate  an  extragalactic  nature  is  image  morphology: 
well-resolved  objects  (as  reported  by  the  photometric  pipeline) 
in  SDSS  are  mostly  galaxies.  However,  these  tests  are  redundant 
in  many  cases,  since  QFOs  in  which  the  host  galaxy  component 
is  strong  enough  to  be  detected  in  the  image  also  generally  show 
noticeable  host  galaxy  features  in  the  spectra.  Furthermore,  these 
tests  work  most  poorly  for  the  most  featureless  objects. 

A  superior  test  is  a  measurement  of  proper  motion,  as  in 
Londish  et  al.  (2002).  Galactic  objects  may  have  measurable 
proper  motions,  depending  on  kinematic  properties,  distance  from 
the  observer,  astrometric  precision,  and  time  baseline.  Since  DC 
WDs  have  low  luminosities,  those  that  are  bright  enough  for 
SDSS  spectroscopy  should  be  relatively  nearby,  within  a  few 
hundred  pc.  For  reference,  a  transverse  velocity  of  20  km  s_1  at 
a  distance  of  200  pc  corresponds  to  a  proper  motion  of  2 1  mas 
vr  1 .  On  the  other  hand,  the  proper  motions  of  extragalactic 
objects  such  as  BL  Lac  objects  should  be  consistent  with  zero. 

As  of  DR2,  proper  motions  for  SDSS  objects  are  determined  by 
matching  with  the  United  States  Naval  Observatory  B1.0  catalog 
(USNO-B;  Monet  et  al.  2003).  In  the  area  covered  by  the  SDSS, 
the  USNO-B  catalog  is  mainly  derived  from  POSS-I  and  POSS-II 
plates,  with  time  baselines  in  the  range  20-5 1  yr.  More  accurate 
proper  motions  can  be  realized  by  using  SDSS  positions  to  recal¬ 
ibrate  USNO-B  astrometry;  the  maximum  time  baselines  ( POSS-1 
to  SDSS)  lie  in  the  range  33-53  yr  for  objects  discussed  in  this 
paper.  This  approach  has  been  implemented  by  Munn  et  al.  (2004) 
and  Gould  &  Kollmeier  (2004).  We  use  the  Munn  et  al.  (2004) 
catalog,  extended  to  include  SDSS  objects  beyond  DR1;  all  but 
four  of  the  580  QFOs  appear  in  the  extended  catalog  at  the  time  of 
writing. 

3.2.1.  Proper-Motion  Assessment 

Some  SDSS  objects  do  not  have  counterparts  in  USNO-B 
because  of  the  difference  in  sensitivity  and  other  factors,  and 
mismatches  occur  in  the  proper-motion  catalog  for  this  and  other 
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reasons.  Therefore,  it  is  necessary  to  evaluate,  based  on  available 
information,  whether  an  individual  proper-motion  measurement 
is  or  is  not  reliable.  We  consider  proper  motions  only  for  objects 
that  are  spatially  unresolved  in  the  SDSS  images;  the  Munn  et  al. 
(2004)  approach  uses  galaxies  (resolved  objects)  to  recalibrate 
USNO-B  astrometry,  and  thus  the  proper  motions  of  resolved 
objects  are  (almost  always)  consistent  with  zero  by  construction. 
We  adopt  the  following  set  of  reliability  criteria:  ( 1 )  we  require  a 
one-to-one  match  between  SDSS  and  USNO-B;  (2)  we  require 
that  the  object  be  detected  on  at  least  four  of  the  five  USNO-B 
plates,  including  both  POSS-I  plates;  (3)  we  require  that  the  rms 
residuals  for  the  proper-motion  fit  be  less  than  525  mas  yr~ 1  in 
both  right  ascension  and  declination;  (4)  we  require  that  there  be 
no  neighbor  with  g  <  22  within  7".  Requirements  (1)  and  (2) 
ensure  that  the  match  between  SDSS  and  USNO-B  is  likely  to  be 
correct.  A  larger  residual  than  required  in  (3)  indicates  that  the  fit 
to  the  proper  motion  was  poor,  in  that  at  least  one  detection  lay 
well  off  a  linear  fit;  525  mas  yr~ 1  is  an  empirical  cut  (more  lenient 
than  that  employed  by  Munn  et  al.  2004  but  the  same  as  used  by 
Kilic  et  al.  2005)  that  includes  the  full  distribution  of  residuals 
seen  for  good  fits.  Requirement  (4)  is  designed  to  avoid  errors 
introduced  by  blending  on  the  Schmidt  plates  that  make  up  the 
USNO-B  catalog  and  is  discussed  in  detail  by  Kilic  et  al.  (2005). 
According  to  these  criteria,  334  of  580  QFOs  have  reliably  mea¬ 
sured  proper  motions. 

In  order  to  confirm  the  utility  of  these  “reliable”  proper 
motions  to  distinguish  between  Galactic  and  extragalactic  ob¬ 
jects,  we  compiled  a  sample  of 3202  spectroscopically  confirmed 
SDSS  quasars  with  secure  redshifis  (and  therefore  no  doubt  as 
to  their  extragalactic  nature).  Of  these,  2092  have  reliably  mea¬ 
sured  proper  motions.  Figure  4  shows  the  distribution  of  reli¬ 
able  quasar  proper  motions;  1987/2092  =  95.0%  have  measured 
proper  motions  g  <  1 1  mas  yr  1 .  This  is  almost  identical  to  the 
Munn  et  al.  (2004)  result,  despite  the  use  of  somewhat  different 
reliability  criteria.  Figure  5  shows  reliable  proper  motions  as  a 
function  of  r-band  magnitude  for  comparison  quasars  and  QFOs. 
It  is  evident  that  a  proper-motion  cut  of  1 1  mas  yr~ 1  effectively 
separates  out  a  large  number  of  QFOs  with  unambiguously  sig¬ 
nificant  proper  motions,  even  toward  the  faint  end.  Therefore,  in 
what  follows  we  shall  consider  proper  motions  g  >  1 1  mas  yr~ 1 
to  be  significant  and  proper  motions  g  <  1 1  mas  yr-1  to  be  in¬ 
significant,  unless  otherwise  stated. 

3.2.2.  Proper  Motions  and  Optical  Colors 

Figure  6  shows  SDSS  color-color  and  color-magnitude  dia¬ 
grams  for  QFOs,  symbolically  divided  according  to  proper  mo¬ 
tion  at  1 1  mas  yr~ 1 .  We  have  not  corrected  for  Galactic  extinction/ 
reddening,  since  many  of  the  QFOs  are  expected  to  be  DC  WDs; 
these  objects  have  low  luminosities  (absolute  magnitudes  Mb  ~ 
14;  McCook  &  Sion  1999),  and  thus  should  be  relatively  nearby 
and  may  not  be  significantly  reddened.  (In  Fig.  6  and  elsewhere, 
we  use  point-spread  function  [PSF]  magnitudes.) 

In  addition  to  a  handful  of  QFOs  scattered  around  in  color 
space,  two  populations  of  QFOs  are  evident  in  Figure  6.  The  ob¬ 
jects  with  significant  proper  motions  are,  as  a  group,  blue  in 
g  —  r,  r  —  i,  and  i  —  z,  and  have  colors  that  approximately  follow 
the  blackbody  color  locus  for  temperatures  in  the  range  7000- 
12000  K.  There  is  also  a  large  group  of  objects  without  signifi¬ 
cant  proper  motions  that  approximately  follow  the  power-law 
color  locus  and  extend  toward  the  red;  these  objects  are  mostly 
concentrated  in  the  range  1  <  aopt  <  2.  (We  define  all  slope 
parameters  a  in  the  sense  fv  <x  v~a.)  These  two  groups  of  objects 
are  remarkably  distinct  from  one  another.  The  region  defined  by 
(g—  r  <  0.35;  r  —  i  <  0.13) — the  “blue  gri”  region — as  de- 


Fig.  4. — Distribution  of  measured  proper  motions  for  2092  spectroscopically 
confirmed  SDSS  quasars.  Proper  motions  are  obtained  from  the  Munn  et  al.  (2004) 
catalog,  extended  to  include  SDSS  objects  not  in  DR1,  and  are  required  to  satisfy 
reliability  criteria  as  described  in  §  3.2.1.  The  vertical  dashed  line  marks  //  = 
1 1  mas  yr-1;  95.0%  of  the  quasars  have  measured  proper  motions  below  this  level. 

lineated  in  the  upper  right  panel  of  Figure  6,  contains  the  vast 
majority  (191/197)  of  objects  with  significant  proper  motions, 
while  it  contains  only  23/137  of  the  objects  without  significant 
proper  motions. 

It  is  plausible  that  the  23  blue  gri  QFOs  with  reliable  proper 
motions  g  <  1 1  mas  yr~ 1  simply  represent  the  low  proper- 
motion  extension  of  the  population  of  QFOs  with  significant 
proper  motions.  The  fraction  (6/120  =  5%)  of  QFOs  outside  the 
blue  gri  region  that  have  reliable  proper  motions  g  >11  mas  yr  1 , 
on  the  other  hand,  is  consistent  with  the  false  positive  rate  for 
quasars  reported  in  the  previous  section.  The  bottom  line  is  that, 
for  objects  with  reliably  measured  proper  motions,  the  blue  gri 
region  contains  almost  exclusively  stars,  whereas  nearly  all  the 
objects  outside  the  blue  gri  region  have  proper  motions  consistent 
with  an  extragalactic  nature.  We  examine  this  proposition  more 
carefully  in  the  following  sections. 

Of  the  246  QFOs  without  reliable  proper  motion  measure¬ 
ments  (including  all  QFOs  with  spatially  resolved  morpholo¬ 
gies),  133  lie  within  the  blue  gri  region  and  113  fall  outside  of  it. 
Thus,  of  580  QFOs,  a  total  of  347  lie  within  the  blue  gri  region 
and  233  lie  outside  of  it. 

The  anticorrelation  between  g  —  r  color  and  /-magnitude  that 
is  apparent  at  the  blue  end,  in  the  lower  right  panel  of  Figure  6,  is 
probably  due  to  a  combination  of  two  selection  effects.  (1)  Very 
blue  QFOs  are  almost  exclusively  stars;  for  the  brighter  objects, 
weak  absorption  features  can  usually  be  seen  in  the  spectra  be¬ 
cause  of  higher  S/N  and  thus  the  objects  are  rejected.  Hence, 
there  are  no  bright  blue  objects  in  the  sample.  (2)  Many  of  these 
objects  are  targeted  for  spectroscopy  by  the  blue  serendipity 
algorithm;  for  this  algorithm  the  most  stringent  faint  magnitude 
cutoffs  are  in  the  u  and  g  bands.  At  a  constant  u  or  g  magnitude, 
the  bluest  objects  have  the  faintest  /-magnitudes. 

3.2.3.  Proper-Motion  Cut 

Within  the  blue  gri  region,  we  reject  all  191  QFOs  with  reli¬ 
able  g  >  11  mas  yr-1  from  further  consideration  as  BL  Lac 
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Fig.  5. — Reliably  measured  proper  motions  vs.  r-band  magnitude  for  comparison  quasars  and  QFOs.  Florizontal  dashed  lines  mark  the  1 1  mas  yr"  1  level;  for  clarity, 
only  500  of  the  comparison  quasars  (randomly  selected)  are  shown. 
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Fig.  6. — Color-color  and  color-magnitude  diagrams  of  SDSS  QFOs;  no  correction  for  Galactic  reddening  has  been  applied.  Contours  represent  the  stellar  locus  (e.g., 
Finlator  et  al.  2000).  Diagonal  magenta  lines  in  the  color-color  diagrams  show  the  loci  of  objects  with  pure  power-law  spectra,  black  dots  on  these  lines  correspond  to 
a  =  —2,  —1,0, 1,2,  and  3  ( fu  oc  v  a)  from  lower  left  to  upper  right,  and  diagonal  black  lines  show  the  loci  of  blackbody  colors.  Red  plus  signs  show  objects  with 
reliably  measured  proper  motions  fj,  <  1 1  mas  yr-1,  blue  crosses  show  those  with  reliable  proper  motions  //  >  1 1  mas  yr-1,  and  gray  triangles  show  objects  that  lack 
reliable  proper  motion  measurements.  The  apparent  correlation  between  g  —  r  color  and  /-magnitude  is  due  to  a  combination  of  SDSS  target  selection  and  S/N  selection 
effects,  as  discussed  in  §  3.2.  The  box  in  the  top  right  panel  is  the  “blue  gri ”  region. 
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Fig.  7. — Similar  to  Fig.  6,  but  broken  up  schematically  according  to  image  morphology  and  redshift.  Gray  triangles  show  objects  with  pointlike  morphologies  and  no 
identifiable  spectral  features  for  obtaining  redshifts.  Blue  plus  signs  show  resolved  objects,  and  red  crosses  show  objects  with  measured  redshifts  large  enough  to  be 
cosmological;  in  many  cases  the  latter  two  types  of  symbol  overlap,  especially  for  the  reddest  objects  (those  with  strong  host  galaxy  components).  The  localization  of 
objects  with  redshifts  or  resolved  morphologies  is  consistent  with  the  separation  according  to  proper  motion  that  is  apparent  in  Fig.  6. 


candidates.  While  the  fraction  of  QFOs  outside  the  blue  gri  re¬ 
gion  that  have  reliable  g  >  1 1  mas  yr~ 1  is  consistent  with  the 
false  positive  rate  for  comparison  quasars,  a  handful  of  objects 
stand  out.  Of  2092  comparison  quasars  with  reliable  proper  mo¬ 
tions,  the  maximum  measured  proper  motion  is  40.7  mas  yr  1 ; 
three  QFOs  outside  the  blue  gri  region  have  measured  proper 
motions  greater  than  this  value  (44.3, 120.3,  and  151.6masyr_1). 
These  three  objects  fall  just  outside  the  blue  gri  region,  all  having 
g  —  r  <  0.35  but  0.13  <  r  —  i  <  0.16.  We  regard  their  proper 
motions  as  genuine  and  reject  them  as  BL  Lac  candidates.  On 
the  other  hand,  the  three  QFOs  outside  the  blue  gri  region  that 
have  reliable  proper  motions  1 1  <  g  <  40.7  mas  yr-1  are  well 
separated  from  the  blue  gri  objects.  They  have  proper  motions 
(13.0,  14.1,  and  23.2  mas  yr-1)  that  are  fully  consistent  with  the 
tail  of  the  quasar  proper  motion  error  distribution.  We  regard 
these  proper  motions  as  dubious  and  therefore  do  not  reject 
these  objects.  Thus,  our  proper  motion  cut  eliminates  a  total  of 
194  QFOs,  leaving  a  total  of  386  BL  Lac  candidates. 

Since  246  of  the  QFOs  lack  reliable  proper  motion  measure¬ 
ments  (and  since  some  DC  WDs  may  have  proper  motions  below 
the  detection  threshold),  the  proper-motion  cut  alone  cannot 
produce  a  clean  sample  of  extragalactic  objects.  Of  the  124 


Kleinman  et  al.  (2004)  objects  that  qualify  as  QFOs  according 
to  our  definition,  39  survive  the  proper  motion  cut  (cf.  §  5.3); 
these  objects  all  lack  reliable  proper  motions. 

We  checked  the  194  QFOs  rejected  by  the  proper-motion  cut 
for  other  properties  suggestive  of  a  possible  extragalactic  /AGN 
nature.  Referring  ahead  to  §  4,  none  of  these  objects  has  a  RASS 
counterpart  within  1',  a  FIRST  counterpart  within  2",  or  an 
NVSS  counterpart  within  10";  thus,  it  appears  that  none  of  the 
QFOs  rejected  by  the  proper-motion  cut  are  radio  or  X-ray 
sources.  In  addition,  none  of  the  QFOs  rejected  by  the  proper- 
motion  cut  has  a  measurable  redshift. 

3.2.4.  Resolved  Objects  and  Objects  with  Redshifts 

Figure  7  shows  the  same  color-color  diagrams  as  Figure  6,  but 
the  different  symbols  highlight  objects  with  resolved  morphol¬ 
ogies  or  convincing  cosmological  redshifts  (either  from  the 
pipeline  or  determined  manually),  that  is,  groups  of  objects  that 
should  be  exclusively  extragalactic  in  nature.  The  near-total 
absence  of  resolved  objects  or  objects  with  redshifts  in  the  blue 
gri  region  is  rather  striking  but  is  in  fact  unsurprising,  consid¬ 
ering  that  significant  host  galaxy  starlight  contributions  (which 
enable  most  of  our  redshift  measurements  and  are  responsible 
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for  the  extended  morphologies)  tend  to  move  objects  toward  the 
red. 

Of  the  347  QFOs  in  the  blue  gri  region,  only  five  have  mea¬ 
surable  redshifts,  and  only  two  are  reported  to  be  resolved  by  the 
SDSS  photometric  pipeline.  The  five  blue  gri  objects  with  red- 
shifts  share  the  following  other  properties :  ( 1 )  redshift  (uncertain 
for  three  out  of  five)  determined  from  weak  broad  emission  lines 
rather  than  host  galaxy  absorption  features;  (2)  unresolved  op¬ 
tical  morphology;  and  (3)  a  proper  motion  that  is  either  unreli¬ 
able  (two  objects)  or  g  <  6  mas  yr~'  (three  objects),  well  below 
our  detection  threshold.  Two  of  these  objects  have  FIRST  coun¬ 
terparts,  and  one  of  these  is  also  a  RASS  source.  These  are  all 
very  interesting  objects,  but  the  (possible)  presence  of  broad 
emission  lines,  although  weak,  already  hints  that  they  may  dif¬ 
fer  in  a  significant  way  from  classical  BL  Lac  objects. 

The  two  blue  gri  objects  with  resolved  morphologies  lack 
FIRST,  NVSS,  and  RASS  counterparts,  and  they  lack  redshifts. 
Visual  inspection  of  the  images  reveals  that  they  are  marginally 
resolved.  These  are  the  only  QFOs  in  the  sample  that  are  spatially 
resolved  but  lack  redshifts;  we  do  not  regard  this  as  convincing 
evidence  of  an  extragalactic  nature. 

4.  CROSS-CORRELATIONS  WITH  X-RAY 
AND  RADIO  CATALOGS 

4. 1 .  ROSAT All  Sky  Survey 

We  cross-correlated  the  full  sample  of  QFOs  with  the 
RASS  bright  and  faint  source  catalogs  (BSC,  Voges  et  al.  1999; 
FSC,  Voges  et  al.  2000).  The  RASS  covers  virtually  the  entire 
sky  in  the  0.1 -2.4  keV  range,  with  a  typical  limiting  sensitivity 
of  10~13  ergs  cm~2  s_1  and  positional  accuracy  of  10"-30"  for 
point  sources.  Using  a  matching  radius  of  1',  we  find  65  matches 
in  the  BSC  and  32  in  the  FSC.  Based  on  a  second  search  using  a 
matching  radius  of  5'  and  applying  a  constant  background  model, 
we  estimate  the  total  expected  number  of  mismatches  within  1'  to 
be  significantly  less  than  1  for  the  BSC  and  approximately  1  for 
the  FSC. 

In  order  to  convert  the  RASS  count  rates  to  approximate  X-ray 
fluxes,  we  use  the  Portable  Interactive  Multi-Mission  Simulator 
(PIMMS;  Mukai  1993).  We  assume  that  the  intrinsic  X-ray 
spectrum  in  the  (observed-frame)  energy  range  0. 1-2.4  keV  is 
well  described  by  a  power  law  with  frequency  index  a  =  1.25 
ex  ^~125),  intermediate  between  typical  values  for  low- 
redshift  LBL  and  HBL  objects  (cf.  Sambruna  1997).  We  modify 
the  power-law  spectrum  according  to  the  Galactic  neutral  hy¬ 
drogen  column  density  /VH  toward  each  object,  obtained  using 
the  FTOOLS12  (Blackburn  1995)  program  NH,  deriving  the 
column  density  values  from  the  maps  of  Stark  et  al.  (1992). 
Adopting  this  absorbed  power  law  as  the  input  spectrum,  we  use 
PIMMS  to  convert  from  observed  ROSAT  PSPC  count  rate  to 
flux  density  at  1  keV  in  the  observed  frame.  For  objects  without 
counterparts  in  either  the  BSC  or  the  FSC,  we  use  the  same 
procedure  to  place  approximate  upper  limits  on  the  1  keV  flux 
density,  assuming  count  rates  just  below  the  FSC  detection 
threshold  (i.e.,  six  counts  divided  by  the  exposure  time  appro¬ 
priate  for  the  object’s  coordinates).13 

4.2.  FIRST  Survey  and  NVSS 

We  also  cross-correlated  the  sample  of  QFOs  against  the 
FIRST  catalog  (Becker  et  al.  1995;  White  et  al.  1997)  and  the 


12  See  http://heasarc.gsfc.nasa.gov/ftools. 

13  An  exposure  map  is  available  at  http://wave.xray.mpe.mpg.de/images/ 
rosat/survey/rass_bsc. 


NVSS  (Condon  et  al.  1998).  The  FIRST  survey  is  designed  to 
cover  approximately  the  SDSS  footprint  area,  with  a  typical 
sensitivity  of  1  mJy  at  1 .4  GHz  for  point  sources,  a  resolution  of 
5",  and  subarcsecond  positional  accuracy.  The  NVSS  covers  the 
whole  sky  north  of  declination  —40°  with  a  resolution  of  45", 
rms  positional  accuracy  better  than  7",  and  a  typical  limiting 
sensitivity  of  2.5  mJy,  also  at  1.4  GHz. 

We  adopted  a  2"  matching  radius  for  FIRST  and  a  10"  radius 
for  NVSS;  optical/radio  matching  is  discussed  in  detail  by  Ivezic 
et  al.  (2002).  We  find  183  and  198  matches  within  these  radii  in 
FIRST  and  NVSS,  respectively.  Of  the  198  NVSS  matches,  32 
are  not  covered  by  FIRST;  the  remaining  166  are  detected  by 
FIRST.  Applying  the  same  type  of  simple  background  analysis 
described  in  §  4.1,  the  total  expected  number  of  mismatches 
using  these  matching  radii  is  near  zero  for  FIRST  and  approxi¬ 
mately  1  for  NVSS.  Our  choice  of  a  10"  matching  radius  for 
NVSS  is  somewhat  conservative;  among  our  QFOs,  there  are 
two  objects  for  which  bright  FIRST  counterparts  (f,  >  50  mJy) 
exist  within  2"  of  the  SDSS  position,  but  the  nearest  NVSS 
counterparts  are  separated  by  >10"  from  the  SDSS  position 
(12"2  and  25"4  separations)  because  of  overlap  with  other  nearby 
bright  radio  sources.  However,  the  potential  for  confusion  be¬ 
tween  SDSS  and  NVSS  sources  increases  rapidly  for  matching 
radii  >10"  and  using  a  matching  radius  of  20"  or  30"  would 
result  in  numerous  false  matches. 

In  order  to  identify  any  likely  FIRST  counterparts  with  optical- 
to-radio  separations  greater  than  the  matching  radius  of  2"  (such 
as  wide-separation  double-lobed  sources),  we  visually  examined 
FIRST  cutout  images  for  all  objects  for  which  FIRST  data  were 
available  at  the  time  of  writing.  Only  one  additional  likely  FIRST 
match  was  discovered.  The  FIRST  source  has  a  somewhat  un¬ 
usual  elongated  radio  morphology  that  approximately  overlaps 
the  SDSS  position  (J231952. 82— 011626.8)  with  a  centroid  sep¬ 
aration  of  3"2;  SDSS  J23 19— 01 16  (which  is  extended  and  shows 
strong  host  galaxy  features)  is  the  only  apparent  optical  source 
in  a  SDSS  mosaic  image  within  5"  of  the  FIRST  centroid  posi¬ 
tion.  There  is  also  a  RASS  BSC  source  with  a  separation  of 
4"  (well  within  the  positional  error)  from  SDSS  J2319— 0116. 
We  consider  SDSS  J2319— 0116  to  be  the  likely  source  of  both 
the  X-ray  and  radio  emission. 

With  very  few  exceptions,  such  as  SDSS  J2319— 01 16,  FIRST 
counterparts  of  SDSS  QFOs  clearly  display  compact  core- 
dominated  morphologies  to  the  limit  of  the  resolution  of  the 
FIRST  survey.  There  are  several  instances  of  objects  showing 
extended  structures  such  as  halos  or  weak  one-sided  jets,  but  no 
apparent  lobe-dominated  cases. 

Where  both  FIRST  and  NVSS  data  exist,  we  use  FIRST  to 
take  advantage  of  its  significantly  finer  angular  resolution.  For  all 
objects  observed  and  not  detected  by  FIRST,  we  use  the  FIRST 
nondetections  to  set  upper  limits  on  flux  density  at  1 .4  GHz  of 
fV' iim  =  5fv  rms  +  0.25  mJy,  where ^ ^  is  the  noise  level  (typi¬ 
cally  «0.15  mJy)  appropriate  to  the  object’s  coordinates  and 
0.25  mJy  is  a  correction  for  the  CLEAN  bias  (e.g.,  Becker  et  al. 
1995).  For  objects  lacking  FIRST  data,  we  do  not  use  NVSS 
nondetections  to  set  upper  limits  because  of  the  potential  for 
confusion  as  discussed  previously. 

5.  PROPERTIES  OF  BL  LAC  CANDIDATES 

Table  1  contains  positions,  SDSS  photometric  data,  redshifts, 
and  proper  motions  for  probable  BL  Lac  candidates.  Table  2 
contains  RASS,  FIRST,  and  NVSS  data  and  derived  quantities  for 
probable  BL  Lac  candidates.  Tables  3  and  4  contain  the  analo¬ 
gous  information  for  the  possible  BL  Lac  candidates. 


TABLE  1 

SDSS  Photometric  and  Spectroscopic  Measurements,  Proper  Motions,  and  Previous  Identifications  for  “Probable”  BL  Lac  Candidates 


IAU  Designation 

a) 

u 

(mag) 

(2) 

9 

(mag) 

(3) 

r 

(mag) 

(4) 

i 

(mag) 

(5) 

Z 

(mag) 

(6) 

a; 

(mag) 

(7) 

Optical  Morph. 
(8) 

Target  Code 
(9) 

Redshift 

(10) 

9 

(mas  yr-1) 

(ID 

Previously  Identified? 
(12) 

SDSS  J012227. 38+15 1023.1 . 

20.34  ± 

0.05 

19.90 

± 

0.02 

19.55  ± 

0.03 

19.25 

± 

0.03 

18.93 

± 

0.05 

0.17 

P 

1048580 

7.1 

N 

SDSS  JO  127 16. 3 1—082128.9 . 

20.26  ± 

0.07 

19.70 

± 

0.03 

19.14  ± 

0.03 

18.81 

± 

0.02 

18.55 

± 

0.04 

0.08 

P 

20 

0.362? 

4.1 

Y 

SDSS  JO  12750. 83— 001346.6 . 

21.03  ± 

0.09 

20.66 

± 

0.04 

19.83  ± 

0.02 

19.30 

± 

0.03 

18.95 

± 

0.04 

0.09 

r 

2097154 

0.4376 

6.0 

N 

SDSS  JO  13408. 95+003 102. 5 . 

21.17  ± 

0.10 

20.22 

± 

0.03 

19.90  ± 

0.02 

19.72 

± 

0.03 

19.79 

± 

0.08 

0.07 

p 

2 

154.5? 

N 

SDSS  JO  14 125.83—092843.7 . 

18.11  ± 

0.02 

17.59 

± 

0.03 

17.21  ± 

0.02 

16.96 

± 

0.01 

16.68 

± 

0.02 

0.08 

p 

7700 

>0.5? 

5.8 

Y 

SDSS  J020 106. 1 8+003400.2 . 

18.90  ± 

0.02 

18.62 

± 

0.02 

18.25  ± 

0.03 

17.99 

± 

0.02 

17.80 

± 

0.02 

0.07 

r 

5636 

0.2985 

7.6 

Y 

SDSS  J020 137. 66+002535.1* . 

20.39  ± 

0.05 

20.06 

± 

0.02 

19.54  ± 

0.02 

19.39 

± 

0.02 

19.38 

± 

0.09 

0.08 

p 

1048576 

176.1? 

N 

SDSS  J022 048. 46— 084250.4 . 

18.92  ± 

0.02 

18.62 

± 

0.03 

18.27  ± 

0.02 

18.03 

± 

0.02 

17.80 

± 

0.03 

0.06 

p 

1056276 

0.5252? 

7.1 

N 

SDSS  J0238 13.68— 09243 1 .4 . 

20.86  ± 

0.09 

20.27 

± 

0.03 

19.63  ± 

0.02 

19.22 

± 

0.03 

18.85 

± 

0.05 

0.08 

r 

2102788 

0.4188 

3.6 

Y 

SDSS  J024 156. 38+00435 1.6 . 

20.66  ± 

0.08 

20.03 

± 

0.02 

19.61  ± 

0.02 

19.41 

± 

0.03 

19.29 

± 

0.05 

0.09 

p 

3 

0.99 

6.1? 

N 

Notes. — Positions  are  in  J2000.0  coordinates.  Photometric  measurements  are  reported  in  PSF  magnitudes,  not  corrected  for  the  Galactic  extinction.  An  asterisk  following  the  IAU  designation  indicates  that  the  SDSS 
data  using  DR2  reductions  were  not  available  for  the  object;  the  data  reported  are  from  the  DR1  reductions.  The  abbreviation  “p”  in  col.  (8)  indicates  that  the  object  is  pointlike  in  the  SDSS  image;  “r”  indicates  that  it  is 
resolved.  Col.  (9)  reports  the  SDSS  primary  spectroscopic  target  code,  in  which  each  (binary)  bit  represents  a  flag  for  a  different  targeting  algorithm;  see  Stoughton  et  al.  (2002)  for  a  list  of  the  different  targeting 
algorithms  and  their  corresponding  bits.  A  question  mark  following  the  entry  in  cols.  (10)  or  (1 1)  indicates  that  the  redshift  is  uncertain  or  that  the  proper  motion  is  deemed  to  be  unreliable,  respectively.  Col.  (12)  indicates 
whether  the  object  was  previously  reported  in  NED  as  a  BL  Lac  object. 
a  Derived  from  Schlegel  et  al.  (1998). 


2552 


TABLE  2 

Multiwavelength  Data  for  “Probable”  BL  Lac  Candidates 


IAU  Designation 

(i) 

JVH 

( 1020cm-2) 
(2) 

ROSAT  All-Sky  Survey 

FIRST 

NVSS 

^ox 

(14) 

(15) 

Match 

(3) 

Sep. 

(arcsec) 

(4) 

Count  Rate 
(S_1) 

(5) 

Exp. 

(s)  fE  (1  keV) 

(6)  (7) 

3  /xb 

(8) 

Sep. 

(arcsec) 

(9) 

/,  (1-4  GHz) 
(mJy) 

(10) 

./rms 

(mJy) 

(ID 

Sep. 

(arcsec) 

(12) 

f,  (1.4  GHz) 
(mJy) 

(13) 

SDSS  J012227. 38+15 1023.1 . 

4.29 

U 

<0.03 

178.2 

<5.7 

<25.7 

0.8 

38.0  ±  1.2 

>0.90 

0.50 

SDSS  JO  127 16. 3 1—082 128. 9 . 

3.97 

U 

<0.01 

438.8 

<3.2 

<14.4 

0.30 

172.2 

0.20 

0.8 

111.9  ±  3.4 

>1.04 

0.60 

SDSS  JO  12750. 83— 001346.6 . 

3.23 

U 

<0.01 

432.2 

<2.8 

<12.3 

0.30 

5.6 

0.15 

4.2 

7.3  ±  0.5 

>0.96 

0.38 

SDSS  JO  13408. 95+003 102. 5 . 

2.91 

U 

<0.02 

373.9 

<2.6 

<11.3 

<1.1 

0.17 

>0.96 

<0.25 

SDSS  JO  14 125. 83— 092843.7 . 

2.66 

F 

19.5 

0.03  ±  0.01 

440.0 

1.5 

6.3 

0.20 

517.8 

0.15 

0.7 

656.7  ±  19.7 

1.45 

0.54 

SDSS  J020 106. 18+003400.2 . 

2.67 

B 

1.9 

0.34  ±  0.03 

400.0 

17.0 

72.1 

0.40 

13.9 

0.14 

2.0 

12.7  ±  0.6 

0.90 

0.34 

SDSS  J020 137. 66+002535.1 . 

2.64 

U 

<0.01 

401.4 

<2.5 

<10.5 

<1.0 

0.15 

>1.02 

<0.22 

SDSS  J022048. 46— 084250.4 . 

3.05 

B 

9.5 

0.09  ±  0.02 

315.0 

4.9 

21.2 

0.20 

57.3 

0.13 

0.9 

64.9  ±  2.4 

1.10 

0.45 

SDSS  J0238 13. 68— 092431.4 . 

2.93 

F 

15.3 

0.03  ±  0.01 

224.0 

1.6 

6.8 

0.60 

3.0 

0.15 

1.08 

0.31 

SDSS  J024 156. 38+00435 1.6 . 

3.45 

U 

<0.03 

226.6 

<3.9 

<17.1 

<1.0 

0.14 

>0.94 

<0.22 

Notes. — All  frequencies  refer  to  the  observed  frame.  In  col.  (3),  ‘ 

‘B”  indicates 

a  match  in  the  RASS  bright 

source 

catalog,  “F” 

indicates  a 

match  in  the  faint 

source 

catalog,  and 

“U”  indicates  that  the  object  was 

undetected,  i.e.,  there  was  no  match  in  either  the  BSC  or  the  FSC  within  1'.  Cols.  (7)  and  (8)  are  derived  using  slightly  different  assumed  power-law  slopes  (see  §§  5.6  and  5.7);  upper  limits  on  these  quantities  are  obtained 
as  described  in  §  4.1.  Lack  of  data  in  cols.  (10)  and  (11)  indicates  that  FIRST  data  were  not  available  for  the  object.  Col.  (11)  gives  the  FIRST  rms  flux  level  appropriate  to  the  object’s  coordinates;  this  is  used  to  set  the 
upper  limits  in  col.  (10)  as  described  in  §  4.2. 
a  10“ 13  ergs  cm-2  s-1  keV-1. 
b  10-13  ergs  cm-2  s-1. 
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TABLE  3 

SDSS  Photometric  Measurements  and  Proper  Motions  for  “Possible”  BL  Lac  Candidates 


u 

9 

r 

i 

z 

Ar* 

A4 

IAU  Designation 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

Optical  Morph. 

Target  Code 

(mas  yr_l) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

SDSS  JO  15604. 20— 010029.0 . 

19.27  ±  0.03 

18.76  ±  0.03 

18.58  ±  0.02 

18.54  ±  0.02 

18.65  ±  0.03 

0.08 

P 

4 

95.3? 

SDSS  J020145.70+01 1013.0 . 

20.97  ±  0.28 

20.14  ±  0.04 

19.99  ±  0.03 

19.97  ±  0.04 

20.12  ±  0.10 

0.07 

r 

2 

SDSS  J0207 12.37+002443.2 . 

20.28  ±  0.06 

20.15  ±  0.03 

20.23  ±  0.03 

20.37  ±  0.04 

20.51  ±  0.16 

0.08 

p 

1048576 

37.3? 

SDSS  J021 142.81+010446.3 . 

18.46  ±0.02 

18.83  ±  0.02 

19.35  ±  0.02 

19.72  ±  0.03 

19.89  ±  0.09 

0.08 

p 

9568256 

0.0 

SDSS  J02 1925.07+004520.3 . 

20.27  ±  0.05 

20.17  ±  0.03 

20.26  ±  0.03 

20.29  ±  0.04 

20.45  ±  0.12 

0.10 

p 

1048576 

6.4? 

SDSS  J03 1303. 3 1—061600.5 . 

20.25  ±  0.05 

20.23  ±  0.02 

20.31  ±  0.03 

20.41  ±  0.05 

20.23  ±0.18 

0.16 

p 

34603008 

8.6? 

SDSS  J03 1729. 73— 071644.7 . 

20.20  ±  0.05 

20.01  ±  0.03 

20.09  ±  0.03 

20.29  ±  0.04 

20.76  ±  0.25 

0.15 

p 

537919488 

SDSS  J033444. 87— 011253.9 . 

20.01  ±  0.05 

19.95  ±  0.03 

19.98  ±  0.02 

20.20  ±  0.03 

20.41  ±  0.12 

0.40 

p 

1048576 

36.2? 

SDSS  J033556. 39— 002425.5* . 

20.73  ±  0.07 

20.31  ±  0.02 

20.11  ±  0.02 

20.03  ±  0.03 

19.97  ±  0.12 

0.34 

p 

2 

19.0? 

SDSS  J034306. 18—054806.3 . 

18.48  ±0.02 

18.36  ±  0.01 

18.48  ±  0.01 

18.63  ±  0.02 

18.83  ±  0.04 

0.14 

p 

538443776 

13.6? 

Notes. — Similar  to  Table  1,  but  for  “possible”  BL  Lac  candidates  (see  §  5.3).  An  asterisk  following  the  IAU  designation  indicates  that  the  SDSS  data  using  DR2  reductions  were  not  available  for  the  object;  the  data 
reported  are  from  the  DR1  reductions.  The  abbreviation  “p”  in  col.  (8)  indicates  that  the  object  is  pointlike  in  the  SDSS  image;  “r”  indicates  that  it  is  resolved.  Col.  (9)  reports  the  SDSS  primary  spectroscopic  target 
code,  in  which  each  (binary)  bit  represents  a  flag  for  a  different  targeting  algorithm;  see  Stoughton  et  al.  (2002)  for  a  list  of  the  different  targeting  algorithms  and  their  corresponding  bits.  A  question  mark  following  the 
entry  in  col.  (11)  indicates  that  the  proper  motion  is  deemed  to  be  unreliable.  None  of  the  possible  BL  Lac  candidates  have  redshift  measurements  or  previous  BL  Lac  object  identifications  in  NED,  so  these  columns  have 
been  excluded. 

a  Derived  from  Schlegel  et  al.  (1998). 
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TABLE  4 

Multiwavelength  Limits  for  “Possible”  BL  Lac  Candidates 


ROSAT  All-Sky  Survey  FIRST 


Na  Count  Rate  Exp.  f,  (1.4GHz) 

IAU  Designation  (1020cm-2)  (s-1)  (s)  //.■  ( I  keV)a  /xb  (mjy)  (mjy)  onx  aIO 

(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 


SDSS  J015604. 20— 010029.0 .  2.57  <0.01  415.7  <2.4  <10.3  <1.0  0.14  >1.17  <0.15 

SDSS  J020145.70+01 1013.0 .  2.81  <0.01  410.2  <2.6  <11.0  <0.9  0.14  >0.95  <0.25 

SDSS  J020712.37+002443.2 .  2.74  <0.02  384.9  <2.5  <10.8  <1.0  0.15  >0.92  <0.27 

SDSS  J021142. 81+010446.3 .  3.10  <0.02  303.5  <2.8  <11.9  <1.0  0.15  >1.04  <0.21 

SDSS  J02 1925.07+004520.3 .  3.17  <0.04  147.4  <5.7  <24.8  <0.9  0.13  >0.79  <0.27 

SDSS  J031303.31— 061600.5 .  5.33  <0.03  226.6  <5.0  <23.5  ...  ...  >0.81 

SDSS  J031729.73— 071644.7 .  5.16  <0.02  271.6  <4.1  <19.2  ...  ...  >0.87 

SDSS  J033444.87-01 1253.9 .  7.51  <0.01  536.7  <4.5  <21.9  ...  ...  >0.91 

SDSS  J033556.39— 002425.5 .  8.00  <0.01  526.8  <4.7  <22.7  ...  ...  >0.88 

SDSS  J034306. 18— 054806.3 .  5.36  <0.01  409.1  <3.8  <17.9  ...  ...  >1.12 


Notes. — Similar  to  Table  2,  but  for  “possible”  BL  Lac  candidates  (see  §  5.3).  None  of  the  possible  BL  Lac  candidates  are  detected  in  the  RASS,  FIRST,  or  NVSS 
catalogs,  so  RASS  match  codes  and  separations,  FIRST  separations,  and  NVSS  separations  and  flux  densities  have  all  been  excluded. 
a  10-13  ergs  cm-2  s  1  keV-1. 
b  10“13  ergs  cn-T2  s“‘. 


5.1.  Redshift  Distribution 

Figure  8  shows  the  redshift  distribution  for  the  154  BL  Lac 
candidates  with  measurable  redshifts  or  lower  limits  from  inter¬ 
vening  absorption;  by  definition  these  all  fall  into  the  “probable” 
category.  The  median  measured  redshift  is  0.45.  It  is  noteworthy 
that  the  redshift  distribution  of  our  BL  Lac  candidates  extends 
out  substantially  farther  than  those  for  previous  BL  Lac  samples 
(e.g.,  Stickel  et  al.  1991;  Rector  et  al.  2000;  Giomrni  et  al.  2005). 
A  handful  of  objects  extend  to  high  enough  redshift  (z^  2.2; 
z max  =  5.03)  that  Lya  enters  the  observed  wavelength  range; 
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these  objects  invariably  show  at  least  hints  of  weak,  broad 
Lya  emission  lines  (typically  with  observed  EW  >  5  A  but 
REEW  <  5  A),  as  in  Figure  3  b.  Some  of  these  high-redshift  BL 
Lac  candidates,  along  with  similar  objects  discovered  in  high- 
redshift  quasar  searches  (Fan  et  al.  1999;  Anderson  et  al.  2001), 
may  in  fact  represent  the  low-EW  tail  of  the  normal  quasar  ( high- 
luminosity  AGN  )  population  rather  than  true  BL  Lac  objects.  All 
of  the  nine  BL  Lac  candidates  with  z  >  2  lack  RASS  counter¬ 
parts,  and  six  of  them  also  lack  FIRST/NVSS  counterparts;  thus, 
it  is  difficult  at  present  to  confirm  or  reject  this  hypothesis  (but 
see  discussion  in  §  5.6).  This  possibility,  in  addition  to  the  large 
sample  size,  may  explain  the  apparent  discrepancy  in  redshift 
distributions. 

Of  the  BL  Lac  candidates  with  no  redshift  information,  those 
that  are  truly  extragalactic  are  likely  to  populate  the  redshift 
range  0.5  Si  z 5/2.2;  lower  than  z~  0.5,  host  galaxy  features 
should  generally  be  apparent  even  in  moderate  S/N  spectra,  and 
atz  >  2.2  the  Lya  forest  would  enter  the  spectra.  Also  shown  in 
Figure  8  is  the  observed  redshift  distribution  of  the  comparison 
quasar  sample,  which  is  close  to  flat  out  to  z  ~  2.  Bearing  in 
mind  that  a  large  fraction  of  the  BL  Lac  candidates  lacking  red- 
shifts  are  likely  to  be  DC  WDs,  it  would  be  impossible  to  make 
the  BL  Lac  candidate  redshift  distribution  match  the  comparison 
quasar  redshift  distribution  up  to  z  ~  2,  even  with  appropriate 
choices  of  redshift  for  the  BL  Lac  candidates  with  unknown 
redshifts  or  redshift  limits.  This  is  qualitatively  consistent  with  a 
picture  in  which  BL  Lac  objects  are  (on  average)  low-luminosity 
AGNs  and  their  higher  luminosity  counterparts,  the  FSRQs,  do 
not  necessarily  share  the  property  of  having  emission  lines  that 
are  weak  with  respect  to  the  beamed  continuum  radiation  (e.g., 
Urry  &  Padovani  1995). 

5.2.  Optical  Colors  of  Objects  with  X-Ray 
and  Radio  Counterparts 


Fig.  8. — Redshift  distributions  for  BL  Lac  candidates  with  measurable  red¬ 
shifts,  as  well  as  the  sample  of  comparison  quasars  from  §  3.2.1.  The  solid  line 
shows  BL  Lac  candidates  with  emission-line  or  host  galaxy  redshift  determi¬ 
nations;  the  dashed  line  shows  BL  Lac  candidates  with  redshift  determinations 
plus  those  with  redshift  lower  limits  from  intervening  absorption.  The  dotted  line 
shows  the  sample  of  comparison  quasars  (scaled  down  by  a  factor  of  approxi¬ 
mately  20).  The  majority  (232/386)  of  BL  Lac  candidates,  including  86/240 
probable  candidates,  are  not  shown  because  of  the  lack  of  redshift  information. 


Figure  9  shows  the  same  color-color  diagrams  as  Figure  6, 
symbolically  coded  according  to  X-ray  and  radio  detections.  Of 
97  RASS  matches,  only  1  lies  in  the  blue  gri  region  of  color 
space.  Of  1 84  FIRST  matches,  only  six  lie  in  the  blue  gri  region 
of  color  space.  Of  32  objects  with  NVSS  matches  and  no  FIRST 
coverage,  only  1  lies  in  the  blue  gri  region.  For  comparison, 
347/580  =  59.8%  of  all  QFOs  lie  in  the  blue  gri  region  of  color 
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Fig.  9. — Similar  to  Fig.  6,  but  broken  up  schematically  according  to  radio  and  X-ray  detections  rather  than  proper  motion.  Contours  represent  the  comparison  quasar 
sample  (see  §  3.2.1).  Gray  triangles  showQFOs  with  neither  FIRST/NVSS  nor  RASS  counterparts.  Blue  plus  signs  and  red  crosses  show  QFOs  with  RASS  orFIRST/NVSS 
detections,  respectively.  (In  many  cases  they  overlap.)  The  QFOs  with  X-ray  and/or  radio  detections  lie  almost  exclusively  outside  the  region  of  color  space  populated 
by  objects  with  significant  proper  motions  (see  Fig.  6). 


space;  after  the  proper-motion  cut,  1 56/386  =  40.4%  of  BL  Lac 
candidates  still  lie  in  this  region.  By  contrast,  there  are  only  22 
QFOs  outside  the  blue  gri  region  that  do  not  have  FIRST  or 
NVSS  counterparts.  This  simultaneously  emphasizes  that  out¬ 
side  the  blue  gri  region  there  is  not  much  room  for  contamination 
by  DC  WDs  and  reinforces  the  notion  that  true  BL  Lac  objects 
with  such  blue  optical  colors  are  much  rarer  than  “normal”  BL 
Lac  objects. 

5.3.  Subdividing  the  Sample 

Primarily  because  we  do  not  have  reliable  proper  motion 
measurements  for  246  QFOs,  the  sample  of  386  BL  Lac  candi¬ 
dates  remains  significantly  contaminated  by  DC  white  dwarfs. 
The  observed  color  separation  in  Figures  6,  7,  and  9,  between 
objects  with  significant  proper  motions  and  objects  with  low 
proper  motions  or  other  evidence  of  an  extragalactic/AGN  nature 
(radio  or  X-ray  detections  or  redshifts),  argues  strongly  that  most 
true  BL  Lac  objects  lie  outside  the  blue  gri  region,  while  most 
DC  WDs  lie  within  it.  Within  the  blue  gri  region  there  are  a  total 
of  10  QFOs  that  show  evidence  of  a  likely  extragalactic/AGN 
nature;  outside  the  blue  gri  region  there  are  a  total  of  11  QFOs 
that  lack  such  evidence  (plus  an  additional  three  rejected  by  the 
proper  motion  cut).  Indeed,  all  available  data  are  consistent  with 


the  more  ambitious  proposal  put  forward  in  §  3.2.2,  that  nearly 
all  the  QFOs  within  the  blue  gri  region  are  stars  whereas  nearly 
all  the  QFOs  outside  the  blue  gri  region  are  extragalactic. 

Based  on  this  compelling  argument,  and  in  order  to  minimize 
the  impact  of  the  remaining  DC  WD  contamination,  we  subdi¬ 
vide  the  sample  into  “probable”  BL  Lac  candidates  and  “pos¬ 
sible”  BL  Lac  candidates.  Here  “probable”  signifies  a  probable 
extragalactic  nature  (though  not  necessarily  classical  BL  Lac 
object  status),  and  “possible”  signifies  a  likely  stellar  nature.  We 
define  probable  BL  Lac  candidates  to  be  those  objects  that  (1)  lie 
outside  the  blue  gri  region  of  color  space  (i.e.,  g  —  r>  0.35  or 
r  —  i>  0.13)  or  (2)  have  X-ray  or  radio  counterparts  or  mea¬ 
sured  redshifts.  We  define  possible  BL  Lac  candidates  to  be  every¬ 
thing  else,  that  is,  objects  that  lie  within  the  blue  gri  region  and 
have  no  indication  of  extragalactic  nature  (X-ray/radio  coun¬ 
terparts  or  redshifts). 

By  these  definitions,  we  have  240  probable  and  146  possible 
BL  Lac  candidates.  Of  the  39  Kleinman  et  al.  (2004)  objects  that 
qualify  as  BL  Lac  candidates,  38  fall  into  the  possible  category; 
the  remaining  object  also  lies  within  the  blue  gri  region  but  has 
an  apparent  FIRST  counterpart.  Thus,  we  conclude  that  the 
probable  BL  Lac  candidates  are  not  appreciably  contaminated  by 
previously  known  WDs.  For  the  possible  BL  Lac  objects,  we 
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offer  a  word  of  caution:  if  there  exist  quasi-featureless  extraga- 
lactic  objects  (no  spectral  features  detected  with  the  SDSS 
spectral  S/N )  that  are  weak  compared  to  classical  BL  Lac  objects 
in  X-ray  and  radio  (undetected  by  RASS  and  FIRST/NVSS)  and 
have  bluer  than  average  optical  colors  (within  the  blue  gri  re¬ 
gion),  they  would  fall  into  the  possible  category.  Although  we 
argue  that  most  if  not  all  of  the  possible  BL  Lac  candidates  are 
likely  to  be  stars,  a  more  detailed  analysis  of  the  statistics  of  DC 
WD  proper  motions  and  number  counts  would  be  required  in 
order  to  prove  that  this  is  the  case.  Such  an  analysis  is  beyond 
the  scope  of  the  present  work. 

5.4.  BL  Lac  Colors  and  Number  Counts 

Adopting  the  conclusion  that  very  few  or  none  of  the  possible 
BL  Lac  candidates  are  in  fact  BL  Lac  objects,  we  note  that  the 
colors  of  the  true  BL  Lac  objects  in  our  sample  (the  probable 
candidates)  are  0.2-0. 4  mag  redder  on  average  than  the  com¬ 
parison  quasars,  shown  as  contours  in  Figure  9.  Starlight  from 
BL  Lac  object  host  galaxies  can  contribute  to  this  effect  if  the 
(apparent)  AGN/host  galaxy  luminosity  ratio  is  lower  for  BL  Lac 
objects  than  for  quasars,  as  expected  if  BL  Lac  objects  are  on 
average  relatively  low  luminosity  AGNs.  Flowever,  subtracting 
the  host  galaxy  contributions  cannot  even  approximately  ac¬ 
count  for  the  differences  in  the  average  colors,  so  we  conclude 
that  the  intrinsic  power-law  (host  galaxy-subtracted)  spectra  of 
our  true  BL  Lac  objects  are  redder  on  average  than  the  intrinsic 
(thermal /disk  plus  nonthermal)  spectra  of  quasars. 

If  we  assume  that  the  probable  BL  Lac  objects  form  a  statis¬ 
tically  complete  sample  to  /  «  19  (just  brighter  than  the  mag¬ 
nitude  limit  for  low-redshift  SDSS  quasar  spectroscopic  target 
selection)  and  take  into  account  that  most  of  the  probable  can¬ 
didates  (all  of  those  with  /  <  19)  have  g  —  i  >  0.2,  then  the 
sample  is  complete  to  at  least  g  w  19.2.  (This  may  be  a  bad  as¬ 
sumption  because  the  sample  may  not  be  complete  or  because  not 
all  objects  in  the  sample  may  be  true  BL  Lac  objects,  but  fortu¬ 
nately  these  two  effects  are  of  opposite  sign.  )  There  are  98  prob¬ 
able  BL  Lac  candidates  with  g  <  19.2,  selected  from  2860  deg2 
of  sky,  yielding  a  surface  density  of  0.034  deg~2  brighter  than 
this  limit.  If  one  is  willing  to  make  the  additional  assumption 
that  g  ~  B,  this  number  may  be  compared  with  the  number  count 
extrapolations  in  Figure  2  of  Padovani  &  Urry  (1991);  they  are 
mutually  compatible,  albeit  with  a  large  range  of  uncertainty. 
For  comparison,  there  are  ~4-5  quasars  deg~2  to  the  same 
limiting  magnitude  (e.g.,  Boyle  et  al.  2000). 

5.5.  SDSS  Target  Selection  Effects 

The  majority  (349/386)  of  our  BL  Lac  candidates,  including 
205/240  of  the  “probable”  candidates,  have  dereddened  u— 
g  <  0.6  (marked  in  the  upper  left  panel  of  Fig.  9),  i.e.,  they 
satisfy  the  UVX  criterion  of  SDSS  low-redshift  quasar  target 
selection.  Of  386  total  BL  Lac  candidates,  225  are  targeted  by 
low-redshift  quasar  target  selection,  including  179/240  probable 
candidates.  The  remainder  of  the  candidates  either  have  colors 
similar  to  WDs  (and  are  therefore  excluded  from  quasar  target 
selection;  Richards  et  al.  2002)  or  are  simply  too  faint.  The  same 
UVX  criterion  provides  grounds  for  inclusion  in  the  blue  ser¬ 
endipity  targeting  routine  as  well;  although  this  selection  route 
receives  lower  priority  than  the  main  target  selection  routines,  it 
populates  parts  of  color  space  that  might  otherwise  be  ignored 
(such  as  the  WD  exclusion  region)  and  extends  to  a  fainter  limit¬ 
ing  magnitude  than  low-redshift  quasar  targeting.  Blue  seren¬ 
dipity  targeting  accounts  for  an  additional  1 06  BL  Lac  candidates, 
including  10  probable  ones.  The  remaining  55  BL  Lac  candidates 


(  including  5 1  probable  ones)  are  mostly  targeted  because  of  their 
radio  or  X-ray  properties. 

5.6.  Radio-Optical  and  Optical— X-Ray  Slopes 

For  the  purpose  of  calculating  broadband  radio-to-optical 
and  optical-to-X-ray  spectral  slopes  (aro  and  aox)  that  can  be 
readily  compared  with  those  appearing  in  other  works,  we  adopt 
standard  (rest  frame)  reference  frequencies  of  5  GHz  =  6  cm, 
6  x  1014  Hz  =  5000  A,  and  2.42  x  1017  Hz  =  1  keV.  Since  we 
do  not  have  observations  covering  these  rest-frame  frequencies 
in  all  cases  (most  notably  near  5  GHz,  where  we  have  essen¬ 
tially  no  data,  the  nearest  observed  frequency  being  1.4  GHz 
from  F1RST/NVSS),  this  procedure  necessarily  requires  some 
assumptions  and  extrapolations.  We  assume  that  the  radio,  op¬ 
tical,  and  X-ray  spectra  are  locally  power  laws  ff  cx  v~a  with 
exponents  ar  =  —0.27  (the  average  value  for  the  1  Jy  sample; 
Shekel  et  al.  1991),  a0  =  1.5  (typical  of  the  probable  BL  Lac 
candidates),  and  ax  =  1-25  (see  §  4.1).  To  objects  that  lack 
redshift  information,  we  assign  the  median  measured  redshift 
(0.45)  if  they  fall  into  the  probable  category,  or  zero  if  they  fall 
into  the  possible  category;  uncertain  redshifts  and  redshift  lower 
limits  are  treated  as  secure  redshift  measurements  for  this  calcu¬ 
lation.  We  normalize  the  assumed  radio  and  X-ray  spectra  using 
the  1.4  GHz  and  1  keV  (observed  frame)  data  points.  After 
correcting  the  optical  magnitudes  for  the  Galactic  extinction 
derived  from  the  Schlegel  et  al.  (1998)  map,  we  normalize  the 
assumed  optical  spectrum  using  the  nearest  in  wavelength  to 
rest-frame  5000  A  of  the  g,  r,  i,  and  z  magnitudes  (depending  on 
redshift). 

Figure  10  shows  an  X-ray-optical-radio  color-color  diagram 
(o:0X ,  aro)  for  probable  BL  Lac  candidates.  There  are  un¬ 
certainties  in  the  locations  of  individual  objects  on  this  diagram, 
mainly  due  to  (1)  nonsimultaneity  of  observations  at  different 
frequencies,  which  is  a  concern  because  of  BL  Lac  variability; 
(2)  uncertainties  introduced  in  the  procedure  of  converting 
RASS  count  rates  to  flux;  and  (3)  nondetection  of  many  objects 
at  X-ray  or  radio  frequencies.  Moreover,  since  we  have  assumed 
uniform  X-ray,  optical,  and  radio  power-law  slopes  to  extrapo¬ 
late  observed  data  points  to  rest-frame  reference  frequencies,  the 
true  scatter  in  the  plot  is  probably  underrepresented.  We  note  that 
Figure  10  qualitatively  differs  from  the  analogous  Figure  3  of 
Collinge  et  al.  (2004),  because  of  the  more  sophisticated  proce¬ 
dure  for  calculating  X-ray  fluxes  and  broadband  spectral  slopes; 
this  simply  serves  to  emphasize  the  potential  for  confusion. 

It  is  nevertheless  apparent  that  the  sample  follows  a  more  or 
less  continuous  distribution  of  properties  from  high-energy 
peaked  (HBL)  toward  low-energy  peaked  (LBL)  SEDs,  with  no 
shortage  of  intermediate-type  objects.  This  behavior  is  charac¬ 
teristic  of  the  blazar  class  (Perlman  et  al.  1998,  2001),  despite 
previous,  and  by  now  historical,  indications  of  a  gap  between 
HBL-  and  LBL-type  objects  (e.g.,  Padovani  &  Giommi  1995), 
which  resulted  from  the  very  different  selection  biases  imposed 
on  radio-  versus  X-ray-selected  samples.  Perhaps  the  most  note¬ 
worthy  feature  of  Figure  10  is  that  all  the  objects  with  X-ray  or 
radio  detections  fit  so  comfortably  within  the  range  of  properties 
of  previously  known  BL  Lac  objects  (e.g..  Fig.  1  ofPerlman  et  al. 
2001),  despite  the  large  sample  size  and  relatively  novel  selec¬ 
tion  biases.  If  there  are  any  BL  Lac  candidates  in  our  sample  with 
exceptionally  weak  radio  or  X-ray  emission,  then  they  must  ap¬ 
pear  in  Figure  10  only  as  limits  or  fall  into  the  possible  category 
(see  discussion  in  §  5.3).  Among  classical  BL  Lac  objects,  the 
lowest  known  values  of  aro  are  in  the  neighborhood  of  0. 1  (cf.  Fan 
et  al.  1999);  four  of  our  candidates  have  upper  limits  on  am  below 
this  level.  Because  our  derived  am  values  are  somewhat  uncertain, 
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Fig.  10. — X-ray-optical-radio  color-color  diagram  for  probable  BL  Lac  candidates.  The  parameters  aox  and  aro  are  defined  in  the  rest  frame  with  the  sign  conven¬ 
tion  fv  oc  v~a,  using  reference  points  at  5  GHz,  5000  A,  and  1  keV  as  described  in  the  text  in  §  5.6.  The  dashed  line  approximately  encloses  the  population  of  known 
BL  Lac  objects  (e.g.,  Fig.  1  of  Perlman  et  al.  2001).  The  slanted  line  shows  the  formal  division  between  HBL-type  and  LBL-type  objects  at  ~  0.75  (e.g.,  Padovani  & 
Giommi  1995).  Red  triangles  show  objects  with  both  radio  (FIRST  =  solid  or  NVSS  =  outline)  and  X-ray  (RASS)  detections;  blue  squares  with  arrows  show  objects 
with  radio  detections  and  X-ray  upper  limits.  The  black  circle  shows  the  lone  object  with  a  RASS  detection  and  a  FIRST  upper  limit.  Pairs  of  arrows  represent  objects 
with  upper  limits  in  both  radio  and  X-rays.  There  are  substantial  uncertainties  in  the  locations  of  objects  on  this  diagram  (see  §  5.6). 


we  separate  out  all  probable  BL  Lac  candidates  with  aro  <0.2 
(eight  objects)  or  no  radio  detections  (another  19  objects).  We 
present  these  27  candidates  in  Table  5;  deeper  radio  obser¬ 
vations  of  these  objects  are  needed  in  order  to  understand  the 
radio-weak  tail  of  the  BL  Lac  object  population. 

In  addition  to  the  broadband  spectral  slope  aro,  another  pa¬ 
rameter  that  has  been  used  to  discriminate  between  radio-loud 
and  radio-quiet  AGNs  is  the  luminosity  density  at  radio  fre¬ 
quencies  (e.g.,  Stocke  et  al.  1992);  these  authors  placed  the  di¬ 
viding  line  at  L„  ~  1033  ergs  s  1  Hz-1  for  z  ~  2  quasars.  This 
number  was  obtained  assuming  an  Hq  =  50  km  s_1  Mpc~'  and 
fltotai  =  flm  =  1  cosmology,  as  opposed  to  the  currently  favored 
“concordance”  Hq  =  70  km  s^1  Mpc^1,  f lm  =  0.3,  and  Q,\  = 
0.7  model.  Conveniently,  these  different  cosmological  models 
yield  nearly  equivalent  (consistent  to  about  2%)  luminosity 
distance  scales  at  z  =  2  (the  concordance  model  gives  a  slightly 
longer  distance),  so  in  practice  no  conversion  is  required.  As¬ 
suming  the  concordance  cosmology,  all  of  the  objects  in  Table  5 
for  which  we  have  radio  constraints  fall  below  the  approximate 
Lv  ~  1033  ergs  s  1  Hz-1  dividing  line.  While  we  do  not  advocate 


this  approach  to  determining  radio-loud/radio-quiet  nature,  this 
analysis  provides  another  useful  benchmark. 

Bearing  in  mind  the  uncertainties  in  our  derived  values  of  am 
and  aro,  HBL  or  intermediate  objects  appear  to  constitute  the 
majority  of  our  sample.  This  may  seem  unsurprising,  consider¬ 
ing  that  the  gap  between  the  optical  and  soft  X-ray  bands  is  less 
than  3  orders  of  magnitude  in  frequency,  while  the  gap  between 
optical  and  radio  is  more  like  5  orders  of  magnitude.  However, 
we  note  that  the  range  of  synchrotron  peak  frequencies  ^peak 
of  BL  Lac  objects  is  thought  to  be  approximately  1013  Hz  < 
^peak  <  1018  Hz  (e.g.,  Anton  &  Browne  2005),  with  LBLs  at  the 
low  end  and  HBLs  at  the  high  end.  Based  on  the  experience  of 
previous  BL  Lac  object  searches — “the  surveys  find  what  they 
can  rather  than  what  exists”  (Padovani  &  Urry  2001) — we  ex¬ 
pected  that  our  sample  would  be  biased  toward  objects  with  ^peak 
close  to  the  optical  (i.e.,  toward  the  low  end  of  the  i/peak  range), 
and  thus  dominated  by  LBLs.  Uncertainties  in  (aox,  aro)  not¬ 
withstanding,  this  is  evidently  not  the  case. 

A  strong  selection  effect  that  may  contribute  to  a  HBL-bias  is  the 
SDSS  quasar  target  selection  UVX  cut  (dereddenedw  g  <  0.6), 


2558 


COLLINGE  ET  AL. 


Vol.  129 


TABLE  5 

Potential  Radio- weak  BL  Lac  Candidates 


IAU  Designation 

(i) 

Optical 

Morphlogy 

(2) 

Redshift 

(3) 

A1 

(mas  yr-1) 

(4) 

^OX 

(5) 

Oro 

(6) 

log10G  QHza 
(7) 

SDSS  J004054.65-09 1526.8 . 

.  P 

5.03 

>1.01 

<0.05 

<32.6 

SDSS  J012155. 87-102037.2 . 

.  r 

0.4695 

8.1 

>1.02 

<0.23 

<30.9 

SDSS  J013408.95+003 102.5 . 

.  p 

154.5? 

>0.93 

<0.29 

<30.9 

SDSS  J020137. 66+002535. 1* . 

.  p 

176.1? 

>0.99 

<0.26 

<30.8 

SDSS  J0241 56. 38+00435 1.6 . 

.  p 

0.99 

6.1? 

>0.90 

<0.22 

<31.5 

SDSS  J024 1 57.37+000944. 1 . 

.  p 

0.7896? 

57.8? 

0.83 

<0.35 

<31.6 

SDSS  J025046. 48— 005449.0 . 

.  p 

10.3? 

>0.74 

<0.28 

<30.8 

SDSS  J0256 12. 47— 001057.8 . 

.  r 

0.6302 

>-0.11 

<0.28 

<31.1 

SDSS  J03 1 712.23—075850.4 . 

.  p 

2.6993 

4.5 

>1.03 

SDSS  J090 133. 43+03 14 12. 5 . 

.  r 

0.4591 

5.0 

>1.04 

<0.21 

<30.8 

SDSS  J  104833.57+620305.0 . 

.  p 

5.0 

>1.07 

<0.27 

<30.8 

SDSS  J1 14153.35+021924.4 . 

.  p 

3.5979 

5.0 

>1.19 

0.11 

32.9 

SDSS  J121221. 56+534128.0 . 

.  p 

3.19 

1.4 

>1.22 

<0.03 

<32.3 

SDSS  J 123743. 09+630 144. 9 . 

.  p 

3.5347 

6.4? 

>1.18 

<0.04 

<32.4 

SDSS  J 124225. 39+6429 19.1 . 

.  r 

0.0424 

5.1 

>1.40 

SDSS  J 1332 19. 65+6227 15. 9 . 

.  p 

3.15 

3.2 

>1.10 

0.19 

32.8 

SDSS  J  142505.61+035336.2 . 

.  p 

2.2476? 

3.6? 

>1.13 

<0.09 

<32.1 

SDSS  J 1508 18. 97+5636 11. 2 . 

.  p 

2.0521? 

3.2 

>0.95 

<0.21 

<32.0 

SDSS  J15 1 1 15.49+563715.4 . 

.  p 

10.6? 

>1.06 

<0.25 

<30.8 

SDSS  J154515. 78+003235.2 . 

.  p 

1.0114? 

4.1 

>1.00 

<0.23 

<31.9 

SDSS  J165806.77+61 1858.9 . 

.  p 

>1.41? 

>1.09 

<0.14 

<31.7 

SDSS  J21 1552. 88+0001 15. 5 . 

.  p 

4.5? 

>0.95 

SDSS  J2 120 19. 13— 075638.4 . 

.  p 

28.7? 

>0.96 

<0.23 

<30.8 

SDSS  J2 13950. 32+104749. 6 . 

.  r 

0.296 

3.2? 

1.02 

SDSS  J224749. 55+134248. 2 . 

.  p 

1.1746? 

5.8 

>1.10 

SDSS  J23 1 000.8 1 —0005 1 6.3 . 

.  p 

>1.68? 

5.0 

>1.04 

<0.13 

<31.9 

SDSS  J232428.43+144324.4 . 

.  p 

1.41 

2.2? 

>0.75 

Note. — Selected  data  from  Tables  1  and  2  for  potential  radio-weak  BL  Lac  candidates  (see  §  5.6). 

a  In  units  of  ergs  s-1  Hz-1,  assuming  Ho  —  70  km  s”1,  —  0.7,  and  =  0.3.  Luminosity  density  has  been  extrapolated  to  rest-frame  5  GHz  as 

described  in  §  5.6.  Uncertain  redshifts  and  redshift  limits  have  been  assumed  to  be  exact,  and  objects  without  redshift  constraints  have  been  placed  at  the 
median  measured  redshift  of  the  sample  (zme d  =  0.45). 


which  most  of  our  BL  Lac  candidates  satisfy,  as  discussed  in 
§  5.5.  Figure  9  shows  that  the  numbers  of  BL  Lac  candidates  drop 
sharply  for  u  —  g  >  0.6.  Thus,  most  of  our  sample  may  be  con¬ 
sidered  UVX-selected.  Although  the  SDSS  u-band  is  centered  at 
t'efif  ~  8.51  xlO14  Hz  (still  toward  the  low  end  of  the  t'peak  dis¬ 
tribution),  the  requirement  imposed  on  UV  continuum  slope  may 
result  in  a  bias  toward  HBL-type  objects. 

5.7.  Comparison  with  Anderson  et  al.  (2003) 

Anderson  et  al.  (2003)  presented  a  sample  of  candidate  BL 
Lac  objects  selected  as  part  of  a  project  to  use  SDSS  imaging  and 
spectroscopy  to  obtain  optical  identifications  of  a  large  sample  of 
RASS  X-ray  sources.  Their  sample  contained38  “probable”  (by 
their  own  definition)  and  seven  additional  “possible”  (by  their 
own  definition)  BL  Lac  objects  selected  from  1400  deg2  of  SDSS 
data.  This  area  of  sky  is  covered  by  the  spectroscopic  parent 
sample  for  our  search.  Of  the  38  Anderson  et  al.  (2003)  probable 
BL  Lac  objects,  our  search  recovers  35  objects;  one  of  the  re¬ 
maining  three  objects  fails  the  S/N  requirement  (see  §3.1.1),  and 
the  other  two  fail  the  polynomial  test  (these  objects  are  classified 
as  galaxies;  see  §  3.1.2).  We  also  recover  two  out  of  the  seven 
Anderson  et  al.  (2003)  possible  BL  Lac  objects.  All  the  Anderson 
et  al.  (2003)  objects  that  we  recover  are  categorized  (according  to 
our  definition)  as  probable  BL  Lac  objects  because  of  their  X-ray 
counterparts. 

Following  Anderson  et  al.  (2003),  we  calculate  the  X-ray-to- 
optical  flux  ratio  fx./fopt,  where  fx  is  the  broadband  absorption- 


corrected  X-ray  flux  in  the  observed  frame  energy  range  0.1- 
2.4  keV and^opt  is  the  broadband  extinction-corrected  optical  flux 
in  the  observed  frame  wavelength  range  4000-9000  A.  As  in 
Anderson  et  al.  (2003),  we  assume  the  X-ray  and  optical  spectra 
to  be  power  laws  fv  ex  u~a  with  ctx  =  1.5  and  aopt  =  0.5.  The 
value  of  ax  assumed  here  differs  by  0.25  from  that  we  employed 
in  §  4.1.  The  choice  of  aopt  =  0.5  is  representative  of  the  sample 
of  QFOs  as  a  whole  (see  Fig.  6),  although  they  span  a  relatively 
large  range  (^3)  in  aopt;  our  probable  BL  Lac  candidates  are 
actually  shifted  toward  steeper  (redder)  slopes  as  discussed  in 
§§  3.2.2  and  5.6.  We  set  the  normalizations  of  the  X-ray  and 
optical  power  laws  using  the  flux  densities  at  1  keV  (cf.  §  4.1) 
and  in  the  g-band  (zxcfr  «  6.36  x  1014  Hz). 

Distributions  of  %Pt  for  several  subsets  of  BL  Lac  can¬ 
didates  are  shown  in  Figure  11.  The  lowermost  histogram  may 
be  compared  with  the  lower  panel  of  Fig.  10  from  Anderson 
et  al.  (2003;  not  shown);  the  Anderson  et  al.  (2003)  distribu¬ 
tion  extends  to  a  slightly  higher  value  of  fx/fopt-  This  proba¬ 
bly  reflects  the  different  selection  biases  at  work;  if  the  eight 
Anderson  et  al.  (2003)  objects  not  included  in  our  sample  are 
fainter  at  optical  wavelengths,  they  will  tend  to  be  biased  to¬ 
ward  X-ray-dominated  objects  for  a  fixed  X-ray  sensitivity.  In 
any  case,  it  is  interesting  that  all  three  distributions  in  Figure  1 1 
and  that  from  Anderson  et  al.  (2003)  peak  at  similar  values  of 
/x//opt,  suggesting  a  broad  similarity  in  the  properties  of  the 
X-ray-selected  Anderson  et  al.  (2003)  sample  and  our  (largely) 
optically  selected  sample.  This  is  consistent  with  the  result  from 
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4.  One  of  the  ways  in  which  BL  Lac  objects  manifest  them¬ 
selves  is  by  their  variability.  SDSS  provides  separate  observa¬ 
tions  of  these  objects  in  the  SDSS  imaging  and  spectroscopy, 
which,  as  Vanden  Berk  et  al.  (2004)  show,  can  be  used  to  explore 
the  statistical  properties  of  the  variability  of  the  population.  Some 
patches  of  the  survey  area  have  multiple  imaging  scans  as  well. 
It  will  be  interesting  to  see  how  our  optically  selected  BL  Lac  can¬ 
didates  compare  in  this  respect  with  previously  known  BL  Lac 
objects. 

5.  Polarization  is  another  defining  characteristic  of  the  BL 
Lac  object  population.  We  have  no  information  about  the  optical 
polarization  properties  of  the  sample,  but  the  objects  are  rea¬ 
sonably  bright,  with  zmed  =  18.95,  so  a  broadband  polarization 
survey  of  a  well-defined  subset  of  the  objects,  especially  those  at 
higher  redshifts,  would  be  straightforward  and  informative. 

6.  Finally,  the  DC  WDs  that  we  have  culled  from  the  sample 
are  interesting  in  their  own  right  (Kleinman  et  al.  2004).  Their 
proper  motions  and  brightnesses  allow  the  calculation  of  crude 
distances  and  thus  a  determination  of  their  luminosity  function 
and  a  comparison  with  the  population  of  cooler  WDs.  This  may 
be  important  for  models  of  white  dwarf  formation  and  cooling. 


Fig.  1 1 . — Histograms  of  log  (fx  /f0 Pt)  for  three  subsets  of  BL  Lac  candidates. 
The  bin  size  is  |  dex,  and  the  middle  and  uppermost  histograms  are  vertically 
offset  for  clarity  by  0.5  and  1.0,  respectively.  The  quantities^  and^pt  are 
calculated  as  described  in  §  5.7.  The  bottom  histogram  may  be  compared  with 
the  bottom  panel  of  Fig.  10  from  Anderson  et  al.  (2003).  The  top  histogram 
consists  of  upper  limits  on  log  (fx/fo pt)  (see  §  4.1)  for  probable  BL  Lac  candi¬ 
dates  undetected  in  the  RASS  BSC  and  FSC. 


the  previous  section  that  the  majority  of  our  objects  appear  to  be 
HBL-type. 

6.  FUTURE  WORK 

The  sample  of  BL  Lac  candidates  presented  in  this  paper  is 
ripe  for  further  analysis  and  future  work.  In  particular: 

1.  None  of  the  BL  Lac  candidates  with  radio  or  X-ray 
counterparts  stand  out  from  the  range  of  properties  of  previously 
known  BL  Lac  objects.  Our  ability  to  constrain  the  relative  abun¬ 
dances  of  anomalously  X-ray  or  radio-quiet  BL  Lac  objects  is 
currently  limited  by  the  sensitivity  of  available  X-ray  and  radio 
observations;  thus,  deeper  observations  at  these  frequencies  are 
needed. 

2.  With  higher  S/N  spectroscopy  of  some  of  the  objects  that 
currently  lack  redshift  measurements,  it  would  be  possible  to  con¬ 
struct  the  luminosity  function  of  BL  Lac  objects;  by  comparing 
this  with  the  luminosity  functions  of  other  radio-loud  AGNs, 
we  could  constrain  models  for  jet  opening  angles  as  a  function 
of  luminosity.  Carrying  out  this  calculation  requires  a  full  un¬ 
derstanding  of  the  selection  function  of  the  sample  imposed  by 
the  SDSS  color  selection  (§  5.5).  For  unreddened  BL  Lac  ob¬ 
jects  atz  <  2,  this  is  reasonably  straightforward,  given  the  high 
completeness  of  the  blue  selection  of  SDSS  quasars  (Richards 
et  al.  2002;  Vanden  Berk  et  al.  2005). 

3.  It  is  also  important  to  quantify  the  incompleteness  of  the 
sample  due  to  our  strict  rejection  of  objects  with  emission  lines. 
There  should  be  BL  Lac  objects  whose  host  galaxies  contain  H  n 
regions,  giving  rise  to  emission  lines  of  high  enough  EW  to  cause 
objects  to  be  rejected  from  the  sample.  There  is  also  some  in¬ 
evitable  incompleteness  at  the  low-luminosity  end  arising  from 
our  requirement  that  the  nonthermal  continuum  be  >50%  of  the 
spectrum  at  either  the  blue  or  red  end,  as  well  as  from  our  re¬ 
quirements  on  S/N. 


7.  CONCLUSIONS 

We  present  a  large  sample  of  candidate  BL  Lacertae  objects 
from  the  spectroscopic  database  of  the  Sloan  Digital  Sky  Survey. 
BL  Lac  objects  are  rare  AGNs  in  which  the  dominant  source  of 
observed  radiation  is  thought  to  be  a  relativistic  jet  of  plasma 
directed  within  a  small  angle  to  the  line  of  sight  and  in  which  the 
broad  emission  lines  typical  of  unobscured  AGNs  are  weak  or 
washed  out.  While  the  latter  trait  is  not  well  understood,  the  up¬ 
shot  is  that  BL  Lac  objects  have  nearly  featureless  optical  spec¬ 
tra;  we  use  this  property  as  the  basis  of  our  selection  criteria.  We 
define  an  object  to  have  a  quasi-featureless  spectrum  if  it  ( 1 )  has 
no  identified  emission  lines  with  rest-frame  equivalent  width 
RFEW  >  5  A,  (2)  shows  no  evidence  of  stellar  absorption  fea¬ 
tures  at  zero  redshift,  (3)  is  not  well  fitted  by  stellar  templates  or 
combinations  of  galaxy  eigenspectra,  and  (4)  has  Ca  n  H  and  K 
break  strength  C  <  0.4  for  |C/(5C|  >  3.  (Certain  extreme  objects 
such  as  some  broad  absorption-line  quasars  formally  fulfill  these 
requirements;  such  “messy”  spectra  are  rejected  as  well.) 

The  primary  astrophysical  contaminants  among  quasi- 
featureless  objects  (QFOs)  are  DC  WDs.  We  have  found  (fol¬ 
lowing  Londish  et  al.  2002)  that  these  can  largely  be  rejected 
based  on  their  high  proper  motions,  although  a  significant  frac¬ 
tion  of  objects  lack  reliable  proper  motion  measurements  and 
hence  a  large  degree  of  contamination  remains  in  our  sample. 
Of  580  originally  selected  QFOs,  the  386  final  BL  Lac  candidates 
are  chosen  to  be  those  objects  with  low  or  unreliable  proper  mo¬ 
tions.  Only  85  of  the  objects  in  the  sample  have  been  previously 
reported  as  BL  Lac  objects  in  the  NASA/1PAC  Extragalactic 
Database  (NED).14  The  sample  is  drawn  from  approximately 
2860  deg1 2 3  of  sky  (~l  of  the  complete  SDSS)  and  has  a  median 
/-magnitude  of  18.95. 

We  have  been  able  to  measure  or  constrain  redshifts  for  1 54/ 
386  of  our  BL  Lac  candidates,  either  from  stellar  absorption  lines 
or,  in  a  few  cases,  by  using  intervening  absorption  systems  to 
place  lower  limits.  The  majority  of  these  objects  (131)  have  red¬ 
shifts  less  than  unity  (the  median  measured  redshift  is  zmed  = 
0.45),  with  a  tail  of  objects  extending  out  to  high  redshift  (zmax  = 
5.03).  These  high-redshift  objects  are  potentially  weak-lined 


14  See  http://nedwww.ipac.caltech.edu. 
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normal  quasars  (as  opposed  to  true  BL  Lac  objects),  like  those 
discussed  by  Fan  et  al.  (1999)  and  Anderson  et  al.  (2001)(butour 
sample  does  not  include  those  specific  objects  because  of  the 
details  of  our  selection  criteria).  The  remaining  true  BL  Lac  ob¬ 
jects  in  the  sample  are  likely  to  have  redshifts  in  the  approximate 
range  0.5  SiziS  2.2. 

We  searched  for  X-ray  and  radio  counterparts  to  QFOs  using 
the  ROSAT All-Sky  Survey  X-ray  source  catalogs  and  the  FIRST 
and  NVSS  radio  catalogs.  On  optical  color-color  diagrams,  fea¬ 
tureless  objects  with  high  proper  motions  tend  to  be  confined 
in  color  space,  generally  following  the  loci  of  blackbody  colors 
as  a  function  of  temperature  (typically  7000-12,000  K).  QFOs 
with  redshifts  or  X-ray  or  radio  counterparts  (i.e.,  evidence  of 
an  extragalactic/AGN  nature)  almost  exclusively  lie  outside  the 
region  of  color  space  containing  high  proper-motion  objects.  In 
fact,  the  optical  colors  of  apparently  extragalactic  QFOs  are  rea¬ 
sonably  well  described  in  terms  of  power-law  spectra  (typically 
1  <  aopt  <  2)  plus  host  galaxy  starlight  rather  than  blackbodies. 
We  conclude  that  most  of  the  156  remaining  BL  Lac  candidates 
with  blackbody-like  colors  are  in  fact  DC  WDs,  either  with  low 
proper  motions  or  too  faint  to  have  reliable  proper  motion  con¬ 
straints.  To  restrict  the  impact  of  this  significant  contamination 
on  our  results,  we  divide  the  sample  into  240  “probable”  (i.e., 
likely  extragalactic)  candidates,  with  colors  unlike  those  of  DC 
WDs  or  other  evidence  of  extragalactic/AGN  nature  (redshift  or 
radio/X-ray  counterparts),  and  146  “possible”  (i.e.,  likely  stel¬ 
lar)  candidates,  with  colors  consistent  with  the  high  proper  mo¬ 
tion  objects  and  no  other  evidence  of  extragalactic/AGN  nature. 
The  probable  BL  Lac  candidates  have  colors  generally  consis¬ 
tent  with  the  ranges  of  known  quasar  colors,  but  on  average  are 
redder  by  ^0.2-0. 4  mag,  plus  a  redder  tail  of  objects  with  more 
significant  host  galaxy  contributions. 

Despite  the  large  sample  size  and  relatively  novel  selection 
biases,  we  find  that  the  multiwavelength  properties  of  our 
probable  BL  Lac  candidates  do  not  systematically  deviate  from 
the  range  of  properties  of  previously  known  BL  Lac  objects. 
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